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ABSTRACT 

We recently described an instability due to the nonlinear coupling of p-modes to (/-modes and, as 
an application, we studied the stability of the tide in coalescing binary neutron stars. Although we 
found that the tide is p-g unstable early in the inspiral and rapidly drives modes to large energies, 
our analysis only accounted for three-mode interactions. Venumadhav et al. showed that four-mode 
interactions must also be accounted for as they enter into the analysis at the same order. They found 
a near-exact cancellation between three- and four-mode interactions and concluded that while the 
tide in binary neutron stars can be p-g unstable, the growth rates are not fast enough to impact the 
gravitational wave signal. Their analysis assumes that the linear tide is incompressible, which is true 
of the static linear tide (the m = 0 harmonic) but not the non-static linear tide (m = ±2). Here we 
account for the compressibility of the linear tide and find that three- and four-mode interactions no 
longer cancel. As a result, we find that the instability can rapidly drive modes to significant energies 
well before the binary merges. We also show that linear damping interferes with the cancellation 
and may further enhance the growth rates. The early onset of the instability (at gravitational wave 
frequencies « 50 Hz) and the large number of rapidly growing modes suggest that the instability could 
impact the gravitational wave signal. Assessing its impact will require an understanding of how the 
instability saturates and is left to future work. 

Subject headings: binaries: close - stars: neutron - stars: oscillations 


1. INTRODUCTION 

Coalescing binaries consisting of two neutron stars or a 
neutron star and a black hole are among the most promis¬ 
ing sources for the new generation of advanced ground- 
based gravitational wave detectors, including Advanced 
LIGO (Harry 2010) in the US, Advanced Virgo (Acca- 
dia et al. 2012) in Italy, and KAGRA (Somiya 2012) 
in Japan. Tidal interactions in such binaries imprint a 
signature on the gravitational waveform that if detected 
would help constrain the highly uncertain neutron star 
equation of state. Quantifying this exciting prospect has 
been the subject of extensive investigation over the past 
three decades. 

The strength of the signature depends on both the am¬ 
plitude of the tidal deformation and its phase lag with 
respect to the line joining each body’s center of mass. Be¬ 
cause the neutron star does not respond instantaneously 
to the changing tidal potential, the gravitational wave- 
induced orbital decay causes a phase lag even in the ab¬ 
sence of viscous dissipation (see e.g., Lai 1994). Indeed, 
numerous studies have shown that if the viscous dissipa¬ 
tion is determined entirely by linear processes, its effect 
on the gravitational waveform is small and the phase 
lag due to orbital decay dominates. These inclucie in¬ 
vestigations of linear dissipation due to fluid viscosity 
(Kochanek 1992; Bildsten & Cutler 1992; Lai 1994) and 
due to the linear excitation of resonant, short wavelength 
internal gravity waves (Reisenegger & Goldreich 1994; 
Lai 1994) and inertial waves (Ho & Lai 1999; Lai & Wu 
2006; Flanagan & Racine 2007). More recent studies 
therefore assume that the phase lag due to viscous dis¬ 
sipation can be neglected during the inspiral (Flanagan 
& Hinderer 2008; Hinderer et al. 2010; Read et al. 2009, 


2013; Damour et al. 2012; Lackey et al. 2012; Lackey & 
Wade 2015). 

As the binary inspirals and the tidal deformation grows 
in amplitude, it becomes susceptible to nonlinear insta¬ 
bilities. These will initially manifest as weakly nonlinear 
wave interactions in which the long length scale tidal per¬ 
turbation excites short wavelength fluid waves within the 
star. These waves remove energy and angular momen¬ 
tum from the tide and thus act as an additional source 
of dissipation not accounted for in linear analyses. Pre¬ 
vious studies argued that such nonlinear effects should 
only become important during the very late stages of the 
inspiral (gravitational wave frequencies /gw ^ 400 Hz) 
because the amplitude of the tidal deformation is too 
small at larger orbital separations. However, the im¬ 
portance of nonlinear wave-tide interactions depends not 
only on the amplitude of the tide but also on the strength 
of the tide’s nonlinear coupling to the internal oscillation 
modes of the star. We must therefore evaluate the non¬ 
linear coupling strengths in order to know when during 
the inspiral the tidal deformation first becomes unsta¬ 
ble and to determine how such an instability affects the 
gravitational wave signal. 

In a previous paper (Weinberg et al. 2012, hereafter 
WAQB), we developed a formalism to study weakly non¬ 
linear wave-tide interactions in close binary systems. 
WAQB focused on the well-known parametric instabil¬ 
ity in which a parent wave (e.g., the tide) resonantly 
excites a pair of short wavelength daughter waves of ap¬ 
proximately half the parent’s frequency (for astrophysi- 
cal applications see also, e.g., Kumar & Goodman 1996; 
Wu & Goldreich 2001; Arras et al. 2003; Weinberg & 
Quataert 2008; Barker & Ogilvie 2010, Essick & Wein- 
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berg 2016). In a follow-up paper, Weinberg et al. (2013; 
hereafter WAB) described a different type of instability 
in which a parent wave nonresonantly excites an acoustic 
wave and a gravity wave (i.e., a p-mode and a p-mode). 
A parent wave that is a p-mode or a tidal perturba¬ 
tion is not resonant with a p-g daughter pair because 
its frequency is much lower than the p-mode’s natural 
frequency. Nonetheless, WAB found that high-order p- 
g daughters with similar wavelength couple so strongly 
that even a highly nonresonant, small amplitude parent 
can excite them. 

The primary application of p-g mode coupling that 
WAB considered was tides in coalescing binary neutron 
stars. They found that the tide is p-g unstable early in 
the inspiral and rapidly drives modes to large energies. 
However, they only considered three-wave interactions. 
Venumadhav et al. (2014; hereafter VZH) showed that 
four-wave interactions enter the analysis at the same or¬ 
der as the three-wave interactions and found, in particu¬ 
lar, that {tide, tide, p-mode, p-mode} interactions cancel 
significantly with {tide, p-mode, p-mode} interactions. 
They concluded that although the non-static tide can be 
p-g unstable, the growth rates are too small to signih- 
cantly influence the inspiral of binary neutron stars. 

The analysis in VZH assumed that the linear tide does 
not compress mass elements in the star, i.e., it is incom¬ 
pressible, and thus V • = 0, where is the dis¬ 

placement field of the linear tide. They relied on this as¬ 
sumption in order to carry out a volume-preserving coor¬ 
dinate transformation that relates the energy of a tidally 
deformed star to that of a radially perturbed spherical 
star. Using this method, they related the four-mode cou¬ 
pling to the three-mode couplings and found that they 
almost perfectly cancel. 

While the static linear tide (the m = 0 harmonic) is 
incompressible, the non-static linear tide (m = ±2) is 
compressible. Its compressibility V • x^^^ increases with 
the ratio of the tidal frequency w = mfl to the buoyancy 
frequency N (where H is the orbital frequency). For an 
inspiraling neutron star w > A^, and the non-static tide 
is highly compressible throughout most of the bandpass 
of gravitational wave detectors such as LIGO, Virgo, and 
KAGRA. 

In this paper, we reanalyze the stability of tidally de¬ 
formed neutron stars to p-g mode coupling. We include 
three- and four-wave interactions and account for the 
compressibility of the non-static tide (and the other finite 
frequency corrections to the static tide). We Hnd that 
these corrections undo the cancellation between three- 
and four-wave interactions and lead to growth rates that 
are faster than the incompressible limit by a factor of 
~ wq/w ^ I, where ojq is the dynamical frequency of the 
star. As a result, there is a large set of unstable modes 
that can potentially reach significant energies before the 
neutron star binary merges. 

We also study how linear damping affects the stability 
of the tide to p-g coupling. Although damping typically 
reduces the growth rates of instabilities, the cancella¬ 
tion between three- and four-wave interactions requires 
a balance between the phases of the p- and (/-mode os¬ 
cillations. We find that the large p-mode damping rates 
counter the stabilizing effects of the cancellation and sig¬ 
nificantly increase the growth rate relative to the inviscid. 


incompressible limit studied in VZH. 

The plan of the paper is as follows. In § 2 we derive 
the nonlinear equations of motion and carry out a sta¬ 
bility analysis in the inviscid, incompressible limit that 
recovers the unstable (but slowly growing) non-static tide 
solution found in VZH. In § 3 we derive the functional 
form of the modified three- and four-wave couplings due 
to the finite frequency corrections (such as compressibil¬ 
ity) . In § 4 we carry out a series of checks in order to test 
the accuracy of our analytic and numerical calculations. 
In § 5 we present the main results of our study, the sta¬ 
bility analysis of tidal p-g coupling with finite frequency 
corrections in coalescing binary neutron stars. In § 6 we 
evaluate how linear damping influences the stability of 
the tide and the p-g growth rate. Finally, in § 7 we sum¬ 
marize and describe the need for a nonlinear saturation 
study. 

2. EQUATIONS OF MOTION 

Let ^(a;, t) be the Lagrangian displacement of the stel¬ 
lar fluid at position x and time t relative to the unper¬ 
turbed spherical background and let the operators 
represent the restoring forces at order i. The equation 
of motion for ^{x,t), including linear forces (/i), three- 
and four-wave nonlinear interactions (/2 and /g, respec¬ 
tively), and tidal forcing (patide) is 

= /iK] + + Patide, (1) 

where p is the background density, 

atide = - VC/ - (I • V) VC/, (2) 

and the £ = 2 tidal potential due to a secondary of mass 
M' in a circular orbit at separation a is 

2 

Uix,t) = -eu;y Y. (3) 

m= —2 

Here e = {M'/M){R/af, M, R, and wq = {GM/Ry/^ 
are the mass, radius, and dynamical frequency of the pri¬ 
mary, H is the orbital frequency, Yim is the spherical har¬ 
monic function, and IU 20 = —\/7r/5, IU 2±2 = v^Stt/IO, 
HQii = 0. Since we restrict the analysis to C' = 2, all 
third and higher derivatives of U in atide vanish. We 
express the perturbation ^ to the spherical background 
as the sum of the linear tide x*'^^ oc e, the second-order 
nonlinear tide x*'^^ oc e^, and a perturbation to the tidal 
flow rj: 

I = + V- (4) 

Keeping terms up to order and linear in r) (because 
we are interested in studying the linear stability of the 
tidal flow to infinitesimal perturbations r/), we have 

P + X^^^ +17)= /i[x^^^] + /i[x^^^] + fiVl] 

+/2[x^^\x^^^] +2/2 [x^^\i7] + 2/2[x^^\'^] 
+3/3[x^^\x^^\i7] 

—pVU — p ^x^^^ • V^ VC/ — p{r] -V) VC/. (5) 
Since the equation of motion of the linear tide is 

dX^'^ =/i[x^'^]-dVC/, (6) 
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and that of the second-order nonlinear tide is 

= fi[x^‘^^] + f2[x^^\x^^^]-p ■ v) VC/, (7) 

the equation of motion for the perturbation rj is 

P»7 = /i W + 2/2 [X^^\ »7] + 2/2 [X^^\ 

+3/3[x^^\x^^^i7]-P(^-V)VC/. (8) 

Expanding in the basis of the star’s linear modes 

^{x,t) = J2vait)^a{x), (9) 

a 

X^^Hx,t)=^Xa\t)^aix), (10) 

a 

X^‘^Hx,t)=^Xa\i)^aix), (11) 

a 

using the fact that the displacements are real, normaliz¬ 
ing the modes according to 

J d^XpCa ■ = EoSab, (12) 

and noting that /i[^a] = ~P^a^a^ obtain the ampli¬ 

tude equation for mode a 

Pa +ujlpa = a;2 ^ [ C/*f, + 2 ^ <bcX^P* 

b \ c 

EiY.<bcdX^c^xT'\ vl ( 13 ) 

C cd / 

where 

Cla6 = -^/d'xp|,-(^,-V)VC/, (14) 

Kabc=^ J ( 15 ) 

Ka6cd = ^y •/a . (16) 

The coefficients are symmetric in all their indices and sat¬ 
isfy the usual selection rules (see, e.g., WAQB) including 
angular momentum conservation in the azimuthal direc¬ 
tion 


we defined the time independent three-mode and four¬ 
mode coefficients 

dX^ab — dJab T 2 ^ ) bCabcXc ^ 

c 

— Uab “1“ 

K4ab = 2 ^ HabcXf''^ + ^ X! biabcdX^c'’X^^ 

c cd 

^(1)ab* 

For analytic purposes, it proves convenient to change 
variables to 

qa = (23) 

Wa 

in which case the amplitude equation becomes 

Qa + [7a - 2imaXi\ Qa + “ {niaXlf' - ima^lXa] qa 

= bOa ^ {Ksab + K^ab) <^bqb- (24) 

b 

This change of variables effectively transforms the ampli¬ 
tude equation from the inertial frame (Equation 20) to a 
frame that is co-rotating with the binary, and thereby in¬ 
troduces Coriolis and centrifugal terms on the left-hand 
side of Equation (24). 


2.1. Characteristic equation of tidal p-g coupling 

In order to determine the stability of the tidal flow to 
an infinitesimal perturbation rj, assume that the ampli¬ 
tude of the latter varies as e*®*. By Equation (24), we 
then have the characteristic equation 

[ -s^ -I- {2maXl -b iXa) S+UjI- {niaXl)'^ - itriaflja] qa 

— X)a ^ ) {dX^ah 4” dX^ab') bObqb: (25) 

b 


or in matrix notation 


-s2 / + {2nL + i'y)s- - iXljL 


-b 


9 = 0 , 


(26) 


where 

(L,7,M(0)) = (ma,7a,w^)^ab, (27) 

= -LdatOb {Ksab + K^ab) . (28) 


m -b rua -b rub = 0, 

(17) 

ma + mb + mc = 0, 

(18) 

nia + mb + nic A md = 0, 

(19) 


for Uabi Habc, and Kabcd, respectively. We can therefore 
write the amplitude equation as 

Va + laPa + Ujlpa=i^l (Ktab + dXlab) 
b 

= ujIJ2 + K4bb) 

b 

( 20 ) 

where the overbar denotes a mode’s complex conjugate, 
the second equality follows because the displacements are 
real, we added a term XaPa to model linear damping, and 


This is similar to Equation (E6) in VZH but with two key 
differences. First, they do not include a linear damping 
term. Second, the nonlinear matrix term in VZH, which 
they write as {0)SA4R{0), assumes that the coupling 
coefficients are the same as those of the static (m = 0) 
tide, modulo a rotation i?(0) by 7r/2 around the y-axis. It 
therefore does not account for Hnite frequency corrections 
to K^ab and K 4 ab (due to, e.g., the compressibility of the 
linear tide). 

Equation (26) is a quadratic eigenvalue problem. To 
reduce it to a standard eigenvalue problem, define 


Ais) 



CB 

I 0 ’ 


(29) 


B = -n^L^-inxL ++ M‘'^\ (30) 

c = 2nL + i'r. (31) 
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Equation (26) can then be expressed as a standard eigen¬ 
value problem A{s)z = 0, i.e., 


C B 
I 0 


2 = SZ, 


(32) 


where z = [sq,q]. In our calculations below, we solve 
this equation numerically using the arpack-|—I- library. 

For analytic study, it will be convenient to rewrite the 
characteristic Equation (25) for the specific case of p-g 
coupling. Letting p subscripts denote p-modes and g sub¬ 
scripts denote p-modes, we have the coupled equations 


- (s - + i-fp^ (s - mp.it) + oj. 


- (s - mg^ fl)^ -b i-ig. (s - mg.il) + 


UJ 




qp, 






^3pi gj^PiQpi + K4g- giiOg^ggi 


(33) 


• (34) 


The coefficients K^p.p. and K^g^g. have magnitude ~ e 
and the coefficient K^p.p. has magnitude ~ (WAQB, 
WAB, VZH). Furthermore, as VZH showed (see their 
Section 2.4), K^p-g. does not enter the stability analysis 
at order e^. Since these coefficients do not alter the sta¬ 
bility at order e^, we do not include them in the analysis. 
Substituting Equation (33) into Equation (34) gives 


- (s - mg.it) + i'jg. (s - mg.it) + LO, 

^4gjgi 

Kspi g^ K^p. giUJp. 


qg, 


gi L 


= ^g, H 

-E 


(s — mp.it) — i^p. (s — mp.ii) — t 


Wg.gg.-(35) 


This characteristic equation determines the stability and 
growth rate of p-g coupling to order e^. 


2.2. Stability analysis in the absence of finite frequency 
corrections and linear damping 

The coupling coefficients K 4 gg and ^ KspgK^pg that 
enter Equation (35) can each individually be of order 
unity or larger for high-order modes. However, if we 
ignore finite frequency corrections, then as VZH showed 

K4g,g. + E Ksp.g.Ksp.g. = (36) 

Pi 


i.e., there is significant cancellation and the residual is 
merely of order ^ 1. Thus, if we ignore finite fre¬ 
quency corrections and linear damping, equation (35) 
gives 


(s-mg^il) - 


OJ] 


gj 


qgj 


= "^9.EE 


Kspigj Ksp.g^ 


(s — mp.il)^ 


gi Pi 


{s — mp.il) —i 


^giqgi^ ( 37 ) 


where we assumed that the cancellation of Equation (36) 
is exact since 0{e^) terms cannot anyways alter the sta¬ 


bility. As in VZH, we can estimate the stability of this 
potentially large system of modes by considering a two 
mode system consisting of a single p-g pair with frequen¬ 
cies ojg il ujp- Since Kspg ~ ecopjojg (WAB, VZH), 
if we define r = s — mpil and a = (jUp — mg)il, the 
characteristic equation of the system becomes 


(r -b a)^ — Wg [r^ — ~ 0. (38) 


There are two stable high frequency solutions near r ~ 
±Wp and two low frequency solutions near 


r ~ —a ± ^ujg — (39) 

If Wg < ejal then the low frequency solutions are unstable 
and the p-g pair grows at a rate 

F ss e |mp — mg \ il. (40) 

By angular momentum conservation (i.e., the coupling 
coefficient selection rules Equations 17-19), only p-g 
pairs that satisfy |mp — mg\ = m couple to harmonic 
m of the £ = 2 linear tide. Therefore, for the static tide 
(m = 0) there is no unstable solution and for the non¬ 
static tide [m = ±2) the growth rate of the instability is 
F « 2eil if ujg < 2eil. This agrees with the results found 
in VZH (see their Appendix E). For a coalescing neutron 
star binary, a growth rate of 2eil is too small to drive p-g 
pairs to large amplitudes before the merger (the modes 
only have enough time to grow by a factor of ^ 10^ in 
energy; see § 5.4). 

In Appendix A we present an argument based on the 
energy principle that shows that the static tide must be 
absolutely stable, consistent with VZH’s result and ex¬ 
tending it to all orders in e. In the remainder of the pa¬ 
per, we therefore focus on the stability of the non-static 
tide. 


3. FINITE FREQUENCY CORRECTIONS 
In their analysis, VZH assume that the shape of the 
non-static linear tide that enters the calculation of 
the three- and four-mode coupling coefficients is exactly 
given by the static (i.e., equilibrium) tide solution 


v-x^'^ = o, 


y(l)=. 

Ar 


u 


(41) 

(42) 


where g is the local gravitational acceleration (as in 
VZH, for simplicity we make the Cowling approxima¬ 
tion throughout our analysis.) However, this solution is 
obtained only in the limit <C N'^ (see, e.g., Terquem 
et al. 1998 and Ogilvie 2014; here uj = mil is the tidal fre¬ 
quency and N is the Brunt-Vaisala buoyancy frequency). 
Otherwise, the linear tide is compressible (V • yf 0). 
In the opposite limit, ^ the linear equations of 
motion are (see § 3.1) 


= -V 



(43) 


and x*'^^ is compressible but irrotational (V x x*'^^ = 0; 
here 6p is the Eulerian pressure perturbation). 

In a neutron star N « O.lojQr/R and for a coalescing 
binary neutron star within LIGO’s bandpass 
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Figure 1. Horizontal displacement of the linear tide at the stellar 
surface (divided by eW 22 li', see Equation (3)) as a function 

of orbital separation a/R. The displacement is found by solving 
the linear inhomogeneous Equation (6) for m = ±2; the peaks 
occur when the driving is resonant with individual £ = 2 modes 

of the star. In Figure 2 we plot V ■ s-^d S't o,/R = 12.5 
and 21, indicated here by the two black squares located between 
resonances. 

Finite frequency corrections to can therefore be sig¬ 
nificant throughout the inspiral and, importantly, not 
only when the orbit sweeps through resonances with in¬ 
dividual modes of the star (i.e., not only when the Lai 
1994 “dynamical tide” is excited^). 

We illustrate this in Figure 1 and Figure 2. Here 
and throughout the paper we use the same neutron star 
model as in WAB; it has a mass M = IAMq and radius 
R ~ 11.7 km and assumes the Skyrme Lyon (SLy4) equa¬ 
tion of state (Chabanat et al. 1998) and a non-buoyant 
crust. In Figure 1 we show the horizontal displacement of 
the linear tide at the stellar surface X^h\R) a function 
of orbital separation a/R. We compute by solving 

the linear inhomogeneous Equation (6) using the method 

described in WAQB. The peaks in Xh\^) ocelli' where 
the linear driving is resonant with an £ = 2 mode of the 
star. In Figure 2 we show the radial profiles of V • x^^^ 
and Xr^'^ for the non-static and static tide at two sepa¬ 
rations not near linear resonances {a/R = 12.5 and 21, 
which for our neutron star model correspond to gravi¬ 


tational wave frequencies /gw — 100 Hz and 45 Hz; see 
black squares in Figure 1). Despite not being near res¬ 
onances, V • x^^^ 3'iid x'r'^ differ significantly from the 
static tide solution due to finite frequency corrections. 

The remainder of this section is organized as follows. 
We first derive (§ 3.1) an analytic estimate of the finite 
frequency corrections to V -X^^^ and X^^^ in order to bet¬ 
ter understand the numerical values of Figure 2. Then, 
over several subsections (§§ 3.2-3.6), we describe our pro¬ 
cedure for calculating the three and four mode coupling 
coefficients. Recall that VZH use a volume preserving 
coordinate transformation in order to map the tidally 
deformed star back into a spherically symmetric configu¬ 
ration and from there calculate the coupling coefficients. 
However, because the linear tide is compressible and thus 
not volume preserving, this procedure cannot be used to 
calculate the finite frequency corrections. Our strategy 
instead is to adopt the more general approach taken by 
WAQB, who calculate the coupling coefficients to order 
e, and to now extend the calculation to order e^. Specifi¬ 
cally, in § 3.2 we consider the finite frequency correction 
to K^pg by making use of the exact expressions for Uab 
and Kabc derived in WAQB. In § 3.3 we derive a sum 
rule for '^KzpgK^pg that allows us to calculate it ac¬ 
curately without doing an explicit numerical sum over 
p-modes. In § 3.4 we describe the four-mode coupling 
coefficient Kabcd that we use to calculate K/^gg. In § 3.5 
we show analytically that the finite frequency corrections 
to the terms that enter at 0{[oJo/u}gY) and ©([wo/wg]^) 
do not undo the cancellation of the zero frequency (i.e., 
incompressible) limit. In § 3.6 we derive expressions for 
the terms that enter 0{\ojQ/ujgY‘). When we numerically 
evaluate these terms in § 5, we find that they do undo 
the cancellation. 

3.1. Estimate of the finite frequency corrections to x^^^ 

In order to obtain an analytic estimate of the finite 
frequency corrections to the non-static linear tide x^^^ 
and its compressibility V • X^^^ consider the equation of 
motion for the linear tide (in the Cowling approximation) 

=-V6p +gSp-pVU, (44) 

where 6p and Sp are the Eulerian pressure and density 
perturbations (here and for the remainder of this section 
we drop the “1” superscript; note that in this section 
p refers to pressure and not the p-mode). Since U ^ 
assume a solution of the form 

x(x,t) = a(x,t)e-*-‘^«, (45) 


^ In Lai (1994), the dynamical tide is defined to be the set of 
resonantly excited modes. In the case of coalescing binary neutron 
stars, these modes have little energy and cannot significantly influ¬ 
ence the inspiral waveform even if they are p-g unstable (WAB). 
However, it is also common to define the dynamical tide as the 
difference in shape between the non-static and static linear tides, 
~ Xe^^ (this is how it is defined by, e.g., Zahn 1970; 
Terquem et al. 1998; Goodman Sz Dickson 1998; WAQB). If the 
finite frequency corrections are large, then x^^^ Xe^^ 
similar and the energy in X^^^ be significant (even between 
resonances). 


where $ = SI. Then 

= o — 2imQct — + imil ) a 


•im4> 


(46) 

Focusing on coalescing binary neutron stars, define the 
orbital decay time due to gravitational wave emission 


ia = TTT = 8.9 
a 


M 


1.2Mq 


-5/3 


f h 


V 100 Hz 


-8/3 


(47) 
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Figure 2. Divergence of the non-static linear tide V • (divided by eW 22 ', left panel) and the radial displacement of the non-static 

linear tide (divided by eW 22 R\ right panel) as a function of radius r at two orbital separations that are not near resonances: the blue 
solid lines are for a/R = 12.5 (/gw — 100 Hz) and the red solid lines are for a/R = 21 (/gw — 45 Hz; see black squares in Figure 1). The 
black dashed lines show the linear static (i.e., equilibrium) tide divergence and radial displacement (divided by €W 2 oR)- 


(Peters & Mathews 1963), where /gw = 
is the gravitational wave frequency and A4 = 

[(MM')^/(M + M')] is the chirp mass {M ~ 1.2Mq 
for M = M' = \AMq). When the driving is not near a 
resonance with a normal mode of the star, the amplitude 
of the non-static tide varies on a timescale ta and (Lai 
1994) 

a ^ a./ta, a ~ (48) 


The duration of a resonance is 5a ~ 0.5i?(a/10i?)“^/"‘ 
(Lai 1994 Equation 3.24, after correcting the typo in the 
exponent for a) and the above scalings should apply when 
outside these narrow resonances (cf. Fig. 1). Since ta » 
17“^, during non-resonant driving 


d'^X _ 2 

^ 


(49) 


perturbation is therefore 


5p = 5p^'^°d -\-uj‘^rpxh 

~ -b w (53) 


where = —pU. Plugging (53) into (52) gives 


Sp ~ + — 

9 


PX 


(eq) 



(54) 


where 5p^'^d) = {U/g)dp/dr. Assuming the perturbations 
are adiabatic, the Lagrangian pressure and density per¬ 
turbations are related by Ap/p = FiAp/p (where Fi is 
the adiabatic index), implying 


5p 5p N'^ 

- = ^ - 1 - Xr 

P Tip g 


(55) 


and we can make the usual separation of variables 

X{x,t) = [Xr{r),rxh{r)^ (50) 

where ui = mXl and we treat the amplitudes Xr and Xh 
as constant in time. 

We can estimate the second order (i.e., w^) finite fre¬ 
quency corrections to with an analysis similar to 
Arras & Socrates (2010). By Equation (44), 


2 5p/p + U 

^ Xh = 

r 


(51) 

2 1 dSp 6p 

dU 

(52) 

^ Xr — -^ 9 -^ 

p dr p 

dr ’ 


In the = 0 limit of the equilibrium tide this gives 
= —U/g and by mass conservation 5p = —V • (px) 
we have 


<1 

II 

o 

(56) 

i.e., 


(eq) _ 

rA2 dr\ J' 

(57) 

We can now solve for the second order correction to Xr 

by plugging (53) into (55) to get 



(58) 


where we expand perturbed scalar quantities, such as 
pressure, as dp{r,9,(j)) = 5p(r)Yim[9,(j)). The cor¬ 
rection to the uj = 0 equilibrium tide Eulerian pressure 


where 


fir) 


1 - 


(eq) 

Xh 

(eq) 

Xr 



_j_d^ 

dr 

Xr 


(59) 
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1 - 


E + 2 
A2 


l + E- 


rN"^ \ 1 (Ping 


) A^ (ilnr^ J 


(here E = £ — ding/dhir). Similarly, we find 

+ I*'”’ 


9 d 

pP dr 


^2 Ar 


(60) 


(61) 


Finally, the second order correction to the divergence is 


V • X 


P 


■ Xr 


dlnp 

dr 


/ OJ ^ 

2 „ r 

y Xr 

ViVy 

Cg A2 


A2 


' w \ 2 g d In p 
^ N / c? d In r 


X 


- (e; + 2)(e; + 1) + 

(eq) 


d? Ing 
din r2 

(62) 


where the last equality is for the specific case of I' = 2 
in the neutron star core, where g ~ ^’KGprjZ. Since 
N ~ 0.1wor/-R and c* ~ l.Swo-R in the core. 


dlnp rN"^ gr gr 

d In r g cl cl 


(63) 


and we find 


V . xci 



— 30 CW 22 



(64) 


We thus see that V • x oc {oj/N)'^ and that for w ~ O.Iwq 
(/ gw ~ 100 Hz), the finite frequency correction to the 
divergence is V • x ~ O.le in the core, consistent with 
the numerical results shown in Figure 2. 


Equations (A55)-(A62) and (A71) in WAQB, 

drpr^'^^gr'V ■ p 


V- T(e + 2) 

Kspg = -^2 


+ ^ ixr - Xh)Pr9h - l^lGpXhPhgh 

+Wp 2^p) XrPhOh + ^pFpXhPrQr 


+KTclV ■ x'^ -PPh 


+TclV -p 


ri + i + 


/SlnTi 

V d\np J g 


rV • xV • g 


+ {xhJ^l - 4xr) V • g - ‘igrV ■ x 
+r ^4^ + V • pXtQt 

+Up [(3Fp + 2E^) Xh + FxXr] Pr9h 

TWp [(3Fg + 2Ep) Xh - 6Txr] Ph9r | + ee>(wg). (66) 


Here g is the (/-mode eigenvector (not gravity), T, Fa, 
and Ga are three-mode angular integrals, and for consis¬ 
tency with the convention used in previous sections we 
write the frequency of the tide as w (rather than uj^). 
For high-order p- and g-modes (Unno et al. 1989; WAB) 


iPr,Ph) 

i.9r,9h) 


Ap 

UJp 

^9 


Cs sin kpr 
cos kpr, - — 

ujpr 


ujg sin kgr cos kgr 


N 


A„ 




, (67) 

( 68 ) 


where Ap^g = yjE^ap^g/pr"^, ap = c^A/ f Cg ^dr, and 
ap = [N/r)/ f Ndlnr. The eO(u!g) corrections consist 
of terms like, e.g.. 


T 

Eg 


drp ( 4,g + r-^ 


’■V • X+ ^^Xr ) Prdr- (69) 
dlnr 


3.2. Finite frequency corrections to K^pg 

As VZH showed and as we found in § 2, the stability of 
p-g coupling depends on the degree of cancellation among 
the terms in 

+ "^2 ^^pgKspg- (65) 

p 

In this section we calculate iXspg accounting for the fi¬ 
nite frequency oj of the linear tide x (here we again drop 
the “1” superscript). WAQB give the full expression for 
Kspg, including finite frequency corrections. The cou¬ 
pling is strongest for high order p-g pairs whose wave¬ 
lengths match (i.e., whose wavenumbers satisfy kp ~ kg), 
which implies LOpfuJo ~ wo/wg. The largest individual 
terms in ATspg are of magnitude ~ e(a;o/wg)^, which we 
write as eO{tJ~‘^). Here we are interested in accounting 
for all terms bigger than or equal to eO(l); since the dy¬ 
namics depend on ^p K^pgKgpg, such terms contribute 
at least at e‘^0{uj~^) and for small enough ujg can lead 
to instability if they are not cancelled by K^gg. From 


A straightforward but lengthy calculation (that makes 
extensive use of integration by parts and the linear equa¬ 
tions of motion) reveals that for high-order p-g coupling, 
Equation (66) reduces to 


^3pg = ^J d'^xpp ■ [{ujI - ^2 xjj^g - iU^Cxg] 


^0 

+eO{uj^ 


where 


(s- 

V)x, 



2 dlnTi 

/SlnTW 1 


dlnp 

V 9\np J g_ 


(70) 

(71) 

V-x^r. (72) 
dr 


For the static tide, a; = 0 and V • x = 0 (i.e., ^^g = 0), 
and our expression for K^pg matches the expression given 
in VZH (see their Equations 54, 62, and 69). The last 
two terms in the integrand of Equation (70) arise due to 
finite frequency corrections to iXspg and come in at order 
eO{l). 
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3.3. A sum rule for ^ K^pgK^pg 

We need to evaluate the finite frequency corrections 
to where the sum is over all p-modes. In 

this section we derive a sum rule for J2p K^pgK^pg that 
reduces the sum to a spatial integral that involves only 
the displacements of the tide and the (/-modes (i.e., the 
p-modes do not explicitly enter the calculation). This al¬ 
lows us to evaluate K^pgK^pg much more accurately 
than if we had to numerically sum over p-modes. For 
simplicity, in this section index a labels the linear tide 
(whose frequency we continue to write as w) and indices 
c and d label the (/-modes. 

To derive the sum rule, we adopt the procedure 
Reisenegger (1994) used to find a sum rule for the linear 
tide overlap integral. First, we expand the vectors '0^^ 
and Cac (Equations 71 and 72) as a sum over modes 

'4’ac ~ Alghc^i,, (73) 

b 

= (74) 

b 

where Mabc and Nabc are the coefficients of the expansion 
which, from the orthogonality of the modes (Equation 
12), are given by 


where the sum is only over p-modes and 


he {p} 

= E 

he {p} 


= E 

he {p} 

, ,2 


1 - 


bOt 


E Ehc7^3h. 

he {p} 

\ 

. .2 1 

M^bc 

_ 

2 abc 

J 

Wfc J 


\ , bO 

1- ^ 

] Mabd - 


/ bO 


%Nabd 

b 


LO 


^abc^abd 2 „ M^^f^^Mabd 


{M:,,Nabd + MabdK^f) 




(w-1) (83) 


(the error is e^O{uj^ ^) because ~ and we 

accounted for all terms bigger than eO{(jJc) in K^bc)- If 
6 is a p-mode then Mabc ~ 0(w“^) and if & is a (/-mode 
then Mabc ~ 0(1). The sum over all modes ^abc^abd 
is therefore dominated by the p-modes, which contribute 
at 0(u;“^); the (/-modes contribute at only 0(1). I.e., 


^ M*abcMabd= Y. M^bcMabd + e^Oil). (84) 
he{p,/,s} help} 


J d^XpCb ■ 4 >ac-. 

d^XpCb ■ Cac- 
Comparing with Equation (70), we have 


Mabc=Y 

P/Q 

Nabc=4- 

rho . 


dA^bc — 


{ujI - Mabc - bjlNabc 


(75) 

(76) 


+ eO(we) (77) 


However, for the sums over all modes M*f^^Mabd/bJb 
and '^M*bcNabd/bj"l, the p-modes and (/-modes both 
contribute at 0(w“^). Since K^pg is large only for 
kb — kc, we show in § 4.4 that we can nonetheless ac¬ 
curately calculate these sums using the approximate re¬ 
lation 


E 

he {p} 


M*abcMabd 


(jJ 


2 

b 


' 2\ ^ ( AI^f^^Mabd 

' be{p,f,g} 


(85) 


for a p-mode ^f,. Let filrj] be the standard linear oscil¬ 
lation operator acting on some arbitrary vector rj. For a 
linear eigenmode Therefore 


d^X^ac ■ fl [^fad]=Y^*^bcMaed f d^X Cb ' /l [C 

be d 


— Eg 'y ) Mabc^abdi 


(78) 


and we find the sum rule 

^ M^abcMabd = -^ ■ fl [^ad] • (79) 

. 0 J 


Using the same procedure, we find the sum rules 


MgbcMabd 

b h 

MabcNabd 


^ J d^Xplpla-lpad 

^ J d^Xpxfla- Cad- 


(80) 

(81) 


Note that these are sums over all modes (p-modes, 
g-modes, and the fundamental (/-) modes). However, 
based on the stability analysis of § 2 (see Equation 35), 
we are interested in the restricted sum 

Klcd + ^ K,b-cK,bd. (82) 

he {p} 


(and similarly for Y.^lbc^abd/bjl), where 


(^h)(r) 


KNcs 

bJcT 


2 


( 86 ) 


which follows from the condition kb — kc- We can there¬ 
fore express the sum over p-modes as a spatial integral 
involving only the tide and the (/-modes 


E dAzbcK^bd 

he {p} 


1 

Eq 


I d^3 


1 - 


2w^ 


^b) 


ll^ac ■ fl bP 


adi 




0 


i^^b) 


ac ' f 1 [Carf] ad * .^1 [C ac \) 


+ e ^{bbc ): 


(87) 


where we used the fact that a — a* since the tidal dis¬ 
placement is real. In Appendix B we show that this result 
can be expressed in the following form (which will prove 
useful for comparing with K 4 cd in § 3.5) 


E dC^bcK^bd 


he {p} 
1 


— y (i x<' Fip 


(aEi) (“Effc) + (oycfi) 


•k s 
^:ks 
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+ (c*a^^) {a%d%) + {d^a%,) 


r 

Tip 


. din g\ 
Wac,r^ad T Wad^r^ac i ^ dlllT j ' 


2lo 


2 d'ipac^r dlpad,^ 
dr dr 


+Wr 


!) I 

f d-tp ac,r dC^ad,r dC^ac^r dlpad,r 


V dr 


dr 


dr dr 


+ e 0{u)^ ), 


( 88 ) 


where the subscript semicolon denotes covariant deriva¬ 
tive and 


Zac{r,0,(j)) = r 


V -xp, 


dxpaCjV 

dr 


(89) 


The first term in Equation (88) is ©(w”^), the next two 
terms are 0{uj~^), and the remaining terms are each 
0{uj~^). In § 4 we compare calculations of 
using both methods — explicitly summing over p-modes 
versus using the sum rule (Equation 87) — and show 
that they are in good agreement. 


Consider a perturbation ^ = a -|- c -|- d, where a is 
the linear tide and c and d are high-order g-modes. The 
ordering of terms in dK 4 /dV is then 

V • adet |(a -I- c-b d)y-| ~(!I(u;“^), 
(V-af 
(a>^) 

V • aV • caLdb + (c o d) ~ ©(w'^), 

V • cdet |(a -b a -b d)L| -b (c O d) ~ 0(u;“^), 

ada^dd‘^p.,jks^O{up‘^), 

(V-a)2(V-c)(V-d)^0(I), 

(V-c)(V-d) a>7~0(I), (94) 

where we assume the strong coupling condition for the 
g-modes kc — kd and thus Wc « ujd- The four-mode 
coupling coefficient, accurate to 0{uj~^), is therefore 


3.4. Four-mode coupling coefficient Katcd 

In order to calculate Kdcd (Equation 22) we need to 
evaluate the four-mode coupling coefficient Katcd- Using 
a Hamiltonian formalism, Van Hoolst (1994) derives an 
expression for Katcd (see Van Hoolst & Smeyers 1993 for 
an alternative derivation). In terms of a displacement 4, 
he shows that the fourth-order coupling density is (Equa¬ 
tion 49 in Van Hoolst with our normalization and in the 
Cowling approximation) 


+3rip(e(,C:') + sripV • 


riP 

6Ao 



(^- 4 ) (V.O" + 3(?^^0 iCa) 

+ 2 ) (V-0'?(,C:' +24V-Idet|?^| 




(90) 


where 


Gi = ri (r? _ 3ri -b 3) -b 3 (Ti - i) 




G 2 = ri (ri- 2 ) + 



s 


(91) 

(92) 


(subscript S indicates derivatives taken at constant en¬ 
tropy). The second equality in the equation for dK^/dV 
follows from the relation (Wu & Goldreich 2001; Schenk 
et al. 2002) 


det kf 




(93) 


^aacd — 

~^ol 


d^xTip 


-bc^ add ’^+ 2V • adet 


-V • cdet 


ijks 

(ci U d '), 


2(«i,c4) (a'dffi^,) + (a>b) (c^ld^) 
(a -b c-b d)b 

(a -b a -b c)!. 


V • d det 


+ (^+ 2 ) ( (V • a)2 cl,d4 +2 (V-a) (V.c)a^d4 


+2 (V • a) (V • d) adjd^ 


+ dO{up^). 


(95) 


This form for Kaacd will be useful for carrying out the 
analytic calculations in § 3.5 at G(a;“^) and G(w“^). In 
Appendix C we give an alternative form for Kaacd that 
is less compact but more suitable for carrying out the 
G(w“^) numerical calculations in § 3.6 and § 5. 


3.5. Finite frequency corrections at O{ojg '^) and 0{u)g ^) 

Having obtained expressions that account for finite fre¬ 
quency corrections to the coupling coefficients to G(a;“^), 
we now evaluate Kdgg + ^ K^pgK^pg and determine 
whether the cancellation (Equation 36) still holds. In this 
section we consider the finite frequency corrections that 
enter at G(w“^) and G(w“^) and show analytically that 
they do not undo the cancellation. In § 3.6 we consider 
the finite frequency corrections that enter at In 

the numerical calculations of § 5 we find that these cor¬ 
rections do undo the cancellation. 

If we again let index a label the linear tide and 
indices c and d label the g-modes, then the G(u;“^) and 
0{uffi) contributions to the coupling coefficient (Equa¬ 
tion 22) 

^4cd — -b cfKaacd (^h) 

come entirely from the G(a;“^) and contributions 

to Kaacd- From the expressions for three-mode coupling 
in WAQB, we know that K^( 2 )ad will at most contribute 
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at 0{uj^ ^). Comparing Kaacd (Equation 95) to our sum 
rule result (Equation 88) we immediately see that the 

terms cancel upon calculating 

^^aacd “t ^ ^3bc^^bd' 

We now show that the 0{uj~^) terms cancel too. The 
remaining terms in K 4 cd + that contribute 

at 0{ijj~^) are, by Equations (88) and (95), 


^biaacd T E K^itcK, 


3bd 


1 

Eq 


(fxr ip 




- V • a det I (a + c + d)]^ \ | + e^O{uj-^) 

• a (flr - a/i) Ch^ (E^ac% - 'S'a.ad.c) 

+db^ {Ea'ld - Sa,ac,d) ] + e^O{u;-^), (97) 


where we used results from Appendix C (specifically. 
Equations C3-C5 and the four-mode angular integrals 
defined in Equations C11-C14). In Appendix C.l we 
show that the angular integrals satisfy the relation 

d^ac,bd ~ ^a,bd,c = ^bc,ad ~ Sa,bc,d- (98) 

The bracketed term in the previous equation is therefore 
proportional to d{chdh)/dr and the expression is seen to 
contribute only at 0{uj~‘^) after integration by parts. 

We have thus shown that the finite frequency correc¬ 
tions that enter Equation (65) at and O{oj~^) do 

not undo the cancellation between three- and four-mode 
coupling. 


3.6. Finite frequency corrections at 0{ujg ^) 

There are a large number of terms that contribute at 
0{ujf^) to K 4 gg + '^KzpgK^pg aud they must be eval¬ 
uated numerically. In Appendix D we describe how 
we carry out the calculation. Eirst, we give numeri¬ 
cally useful expressions for the terms in 3K^(i)^(i)gg + 
EspgKspg. We then consider two methods for calcu¬ 
lating the nonlinear tide which we need in order 

to compute 2K^{2)gg. The first method involves expand¬ 
ing as a sum over modes and solving for the co¬ 
efficients of the expansion. The second method involves 
directly solving the inhomogeneous equation of motion of 
the nonlinear tide (Equation 7). We find that the second 
method is significantly more accurate. We will present 
the results of the ©(wj^) calculation in § 5. 

4. TESTING THE CALCULATION 

In order to test the accuracy of our analytic and numer¬ 
ical calculations, in this section we perform the following 
four checks: 


1. In § 4.1 we evaluate whether our numerical cal¬ 
culations can recover the expected cancellation to 
©(e^) in the case of the static tide. 

2. In § 4.2 we compute 3«:^(i)^(i)gg -f X) KspgKzpg by 
explicitly summing over p-modes and integrating 


the full four-mode coupling coefficient. We then 
compare this with the analytic method described in 
§§ 3.5 and 3.6 in which the cancellation at ©(u;”^) 
and is first carried out by hand and then 

the remaining 0{u}f‘^) contributions computed nu¬ 
merically. 

3. In § 4.3 we check our solution of the nonlinear tide 

by comparing results from our two calcula- 
tional methods (solving for the coefficients of the 
mode expansion versus directly solving the inho¬ 
mogeneous equation). 

4. In § 4.4 we compute the sum over I/w^ terms in 

Equation (83) by directly summing over p-modes 
and compare this result with our approximate an¬ 
alytic expression for the sum (the terms in 

Equation 88). 


4.1. Testing the static tide cancellation 
If we only include the m = 0 component of the linear 
tide in our calculation, then there are no finite fre¬ 
quency corrections and based on the analysis in VZH and 
the argument from energy principles given in Appendix 
A, we expect near-exact cancellation (i.e., cancellation to 
0{e^)-, see Equation 36). We have already shown analyt¬ 
ically that the and 0(ujg^) contributions cancel. 

We now want to test whether in the case of the m = 0 
static tide our numerical calculation of the terms 

cancel as expected. To provide a measure for the degree 
of cancellation, define the residual 

^99 ~ + E/ d^3pg^3pg 

p 

— 2K^(2)gg + 3K^(l)^(l)gg A j2K,pgK,pg (99) 

P 


and the fractional residual agg 


^aa — 


TZa 


2K-y(2) f 


( 100 ) 


For the case of the static tide, we expect TZgg ~ 0{e'^) 
and agg ~ (wg/wo)^ since K^( 2 )gg ~ e‘^{uJo/ujg)'^. 

To carry out this calculation, we only include the m = 
0 linear tide and compute the coupling coefficients us¬ 
ing Equations (D11)-(D18) and the direct solution 
(Equation D30). We find the p-modes (and the p-modes 
in the calculations below) using the Aarhus adiabatic os¬ 
cillation package ADIPLS (Christensen-Dalsgaard 2008). 
At the maximum resolution, our background model and 
eigenfunctions have ~ 5 x 10^ grid points over the star. 

In Figure 3 we show agg for the static tide coupled to 
high-order £ = 4 g-modes with 100 ^ n < 1000. We 
find agg ~ 10“® over the entire range in n, implying that 
the individual terms in TZgg cancel to a part in ^ 10® in 
our calculation. For these modes, ujg « 0.2uJo/n and we 
therefore expect 10“® < agg < 10“®. However, given the 
resolution set by the 5 x 10® background grid points, a 
value of agg ~ 10“® is consistent with near-exact cancel¬ 
lation (we checked that agg increases in proportion to de¬ 
creasing background resolution). We therefore conclude 
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Figure 3. Fractional residual agg as a function of ^-mode radial 
order n for the case when only the static (m = 0) tide is included. 
Each point is for an £ = 4 daughter pair of the neutron star model 
with ni = 712 and azimuthal orders that must satisfy mi -|-m 2 = 0. 



that our calculation successfully recovers the correct re¬ 
sult in the static tide (i.e., incompressible) limit. 

4.2. Testing the calculation of 3K^^gg + J2 ^Spg^Spg 

In order to accurately calculate 

3K^(i)^(i)gg + 'y ^ KspgKspg ( 101 ) 

P 

we carried out a series of analytic steps, which we de¬ 
scribed in detail in § 3. Briefly summarizing these steps, 
we showed in § 3.3 that the sum over p-modes can be 
expressed as an integral that involves only the tide and 
the p-modes. After simplifying that integral (Appendix 
B), we showed in § 3.5 that the two cancel to 0(a;“^) and 
even with finite frequency corrections. We then 
rewrote the contribution in a more numerically 

accurate form (Equations D11-D18). 

As a check of these analytic steps and their numeri¬ 
cal implementation, we compare our results to a more 
direct, though less numerically accurate, calculation of 
Equation (101). We carry out the comparison for the 
full tide, including the non-static components, in order to 
test the calculation with the finite frequency corrections. 
In Figure 4, the blue solid line shows the results given by 
the analytic approach summarized in the preceding para¬ 
graph. The red dashed line shows the results given by 
the explicit sum over p-modes and relying entirely on the 
numerics to handle cancellations between terms (includ¬ 
ing the very large 0{ujg‘^) and O(ujg^) contributions). 
We describe the latter calculation in the next paragraph. 
We see that the red and blue lines nearly overlap over 
the full range in m. While the agreement is very good, it 
is not perfect. Given that the analytic approach passes 
the static tide test to ~ I0“® precision (§ 4.1), the most 
likely explanation for the small but > 10 “® differences 


between the blue and red lines is numerical imprecision 
on the part of the red line, which must rely entirely on 
the numerics to cancel the large O(ujg^) and O(ujg^) 
terms. This test should therefore give us confidence in 
our analytic approach to calculating Equation (101) and 
its numerical implementation. 

The p-modes that contribute to the sum are those for 
which kp ~ kg, which corresponds to uip ~ AgNcg/tOgr. 
For a typical stellar model, including our neutron star 
model, the quantity Ncs/r varies somewhat with radius, 
especially outside the stellar core. As a result, there are 
many p-modes that contribute to the sum. In order to 
minimize the error that arises from summing over many 
modes, rather than use a neutron star model to carry out 
the test, we use a polytropic model and set the adiabatic 
index Ti (r) such that 



where A is a constant. Specifically, when calculating the 
poly trope’s buoyancy 


we set 


N^ = -g 


dlnp 

dr 


1 dlnp 
ri(r) dr 


ri(r) = 7 



(103) 


(104) 


where 7 = 1 - 1 - 1 /ripoiy and ripoiy is the (constant) poly¬ 
tropic index. The results in Figure 4 are for an £ = 4, 
n = 50 p-mode of an ripoiy = 2 polytrope with A = O.Iwq, 
the latter being similar to the value in the core of our 
neutron star model. In Figure 5 we show which 

enters the calculation of ^ K^pgK^pg, as a function of the 
radial order Up of the p-modes included in the sum. The 
magnitude of K-^{i)pg peaks sharply for those p-modes 
that best satisfy the condition kp = kg. For modes that 
are Aup ~ 10 away from the peak, K^Mpg is ~ 10 ® times 
smaller and therefore the p-modes in the tail make only 
a very small contribution to ^ KspgK^pg (if we had not 
set A to be a constant the distribution would be much 
broader and many more p-modes would contribute). 


4.3. Testing the calculation of K^( 2 )gg 

Here we compare our direct solution of the nonlinear 
tide found by solving the inhomogeneous Equation 
(D30) with the mode expansion solution found by solv¬ 
ing Equation (D23). While we expect the former to be 
more numerically accurate, such a check is useful, es¬ 
pecially given the complicated forcing terms (Equations 
D39, D40). The black and magenta curves in Figure 4 
compare K^( 2 )gg calculated using the two different meth¬ 
ods for the same polytropic model and modes as § 4.2. 
We see that they are in good agreement. Given that the 
high-precision static tide cancellation made use of the 
direct solution of X^^\ the small differences for certain 
values of m are more likely due to some imprecision in 
calculating x^^^ using the mode expansion. 

We find in § 5 that the finite frequency corrections 
due to the non-static tide do not cancel and yield agg ~ 
O.OI — O.I (when not near a linear resonance). While 
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Figure 4. Comparison of coupling coefficients calculated using the 
analytic approach and the direct numerical approach. We show 
results for the i = 2, n = 1000 p-mode pairs of the npoiy = 2, 
A = O.IcJq polytrope at afR= 1.23 (corresponding to the /-mode 
resonance). The first three points along the ic-axis correspond to 
the p-mode pairs with mi -\-m 2 = 4, the next four points are those 
with mi -|- m 2 = 2, and the final five points are those with mi -h 
m 2 = 0. The blue solid line shows the “analytic” calculation of 
Equation (101) as given by the sum of the eight terms discussed in 
Appendix D.l. The red dashed line shows the numerical calculation 
of Equation (101) as given by explicitly summing over p-modes. 
The black solid line shows the calculation of 2k ( 2 ) with 
found by directly solving the inhomogeneous Equation (D30). The 
magenta dash-dotted line shows the calculation of ‘^‘^■^{‘ 2 )with 

found by using the mode expansion solution (Equation D23). 
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Figure 5. Three-mode coupling coefficient ^■^ii)pg ^ function 
of the p-mode radial order rip for the p-modes that dominate the 
sum given by the blue line in Figure 4. 


the checks above suggest that the analytic approach is 
in good agreement with the direct calculation involv¬ 
ing sums over modes, the limited accuracy of the latter 
precludes us from using it to calculate such small crag’s. 
However, near resonances the residual agg becomes espe¬ 
cially large and even with the direct calculation it is easy 
to distinguish 2 K^( 2 )gg from 3K^(i)^(i)gg -I- 
and see that the 0{uj~'^) terms do not cancel as a re¬ 
sult of the finite frequency corrections. The calculation 
shown in Figure 4 is at the resonance with the /-mode 
(a/R = 1.23) of the ripoiy = 2 polytrope. We see that 
the residuals calculated using the analytic approach and 
the direct approach are in good agreement, with both 
yielding agg ^ 100 for the m = 2 and m = 4 harmonics 
of the nonlinear tide. 


4.4. Testing the 1/w^ sum rule terms 

As we describe in § 3.3, when calculating ^3pg^3pg 
using the sum rule (Equation 87), we approximate the 
terms as 


E 

be {p} 


MXMabd 


(jJ 


2 

b 


' 2\ E! M^f^^Mahd 

' b(i{pJ,g} 


1 

Eq 


^3^'4’*ac- fli^ad] 

LL JL / o \ 


(105) 


(and similarly for the Wj, ^ terms involving Nate), where 

(loI) = (AgNcs/ujgr)'^. This approximation enters our 
calculation of K^(i)^{i)gg+Y^ K^pgK^pg as term IV of Ap¬ 
pendix D (Equation D14). In order to check the accuracy 
of the approximation, we explicitly solve the sum over p- 
modes using the same polytropic model and g-mode pairs 
as § 4.2. We show the results in Figure 6 for a/R = 1.23 
(near the £ = 2 /-mode resonance) and a/R = 25 (not 
near a resonance). We hnd good agreement between the 
analytic calculation given by term IV (blue curves) and 
the numerical sum over modes (red curves). At large or¬ 
bital separations, term IV makes only a very small con¬ 
tribution to the total coupling K^(i)j^(i)gg + ^^ 3 pgl^ 3 pg 
(which is shown as the black curve). At smaller separa¬ 
tions, the contribution is more significant but still sub¬ 
dominant and much less than 2 K^{ 2 )gg (cf. black curve 
in Fig. 4). As a result, we find that the fractional er¬ 
rors introduced by the approximations in Equation (105) 
are always much smaller than the values of the residual 
agg (~ 0.01 — 0.1 when not near a resonance) and there¬ 
fore conclude that they do not affect the accuracy of our 
calculation. 


5. RESULTS WITH FINITE FREQUENCY CORRECTIONS 

In this section we carry out a stability analysis of p-g 
coupling accounting for the finite frequency corrections 
due to the non-static tide of a coalescing binary neutron 
star. In § 5.1 we calculate the degree of cancellation 
between three- and four-wave coupling as measured by 
the residual TZgg (Equation 99). In §§ 5.2 and 5.3 we 
evaluate the stability and growth rate of p-g mode pairs 
given the computed values of TZgg- We use these results 
in § 5.4 to estimate the number of e-foldings of growth an 
unstable p-g pair can undergo before the binary merges 
(ignoring the nonlinear saturation of the instability). 
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Figure 6. Comparing the two different methods for calculating the 1/^^ coupling terms of K^pgK^pg given by Equation (83). As in 
Figure 4, we show results for the £ = 2, n = 1000 ^-mode pairs of the Upoiy = 2, A = 0.1a;Q polytrope. The left panel is for ajR = 1.23 
(corresponding to the /-mode resonance) and the right panel is for a/R = 25 (which is not near a linear resonance). The first three 
points along the x-axis correspond to the p'-mode pairs with mi -|- m 2 = 4, the next four points are those with mi -|- m 2 = 2, and the 
final five points are those with mi -h m 2 = 0. The blue solid line shows the analytic method for calculating the terms (which 

makes the approximation described in § 3.3 and is given by term IV, Equation D14). The red dashed line shows the numerical method 
for calculating the 1/^^ terms (explicitly summing over p-modes). For comparison, the black solid line shows the analytic method for 
calculating ^ ^3pg^3pg given by the eight term sum of Equations (Dll)—(D18). 


5.1. Residual 7Zgg due to finite frequency corrections 

If we ignore finite frequency corrections then we ex- 
pect a residual TZgg ~ 0{e^). By Equation (100), this 
corresponds to a fractional residual agg ^ {uig/ujop since 
K.x^'tgg ~ ^‘{ujQ/wgY. For illustrativc purposes, in this 
section we show results from our SLy4 neutron star 
model for 5 -mode pairs ( 51 , 52 ) with £ = 2, n = 1000 , 
and |mi + 1712 ] = {0,2,4} (other combinations do not 
couple to the tide). The results we show are represen¬ 
tative of what we find with other high-order, low £, 5 - 
modes. Since this particular 5 -mode has a frequency 
ojg ~ 10 “‘*wo, if finite frequency corrections to the cou¬ 
pling coefficients are not significant we would expect 
agg ~ 10“®. However, based on the results of § 4.1, we 
can only expect our calculation of agg to be accurate at 
the ~ 10 “® level (as determined by the resolution of our 
model and eigenfunctions). Therefore, if the computed 
agg < 10 “® it suggests finite frequency corrections are 
not significant and vice versa. 

In Figure 7 we show agg accounting for finite frequency 
corrections due to the non-static tide. The red curve 
shows agg as a function of orbital separation a/R. The 
blue curve is the same as that of Figure 1 and illustrates 
the variation of with a/R and the location of linear 
resonances. The 5 -mode pairs in the three panels are 
i = 2, n = 1000 , and (mi,m 2 ) = ( 2 ,- 2 ), ( 2 , 0 ) and 
( 2 , 2 ); other possible pair combinations such as ( 1 , 1 ) and 
( 1 ,- 1 ) yield similar results. 

We find that between resonances, 10“^ < agg < 1 for 
10 ^ /gw ^ 10^ Hz.^ This result suggests that the fi- 

^ Although agg can be much larger near resonances, because 
our calculation does not account for the limited duration of each 


nite frequency corrections are significant and undo the 
cancellation of three- and four-wave coupling found in 
the static tide (i.e., incompressible) limit. The vari¬ 
ations in agg with a/R are due to variations in 
and x^^^ with a/R and, as the figures show, correlate 
with the locations of linear resonances. The magnitude 
of agg tends to decrease with increasing a/R, as might 
be expected given that finite frequency corrections, in¬ 
cluding the compressibility V • X^^\ decrease with de¬ 
creasing driving frequency. Note, however, that although 
V 'X^^^ oc {w/wqY oc {R/ap (see Equation 64), the mag¬ 
nitude of agg falls off more slowly than [R/ap. As we 
explain more fully in § 5.3, this is because even at larger 
separations, the finite frequency corrections remain sig¬ 
nificant near the inner turning radius of the tide (where 
u) ~ N{r)). Since the nonlinear coupling is strongest 
near the inner turning radius, the residual approaches 
the static tide limit more slowly than [R/ap. 

5.2. Stability analysis 

We return to the stability analysis of § 2 but now in¬ 
clude the effects of finite frequency corrections. From the 
definition of TZgg (Equation 99) 

— ~ ^3pigi +T^gjgi- (106) 

Pi 

The characteristic equation for p-g coupling (Equation 
35) including finite frequency corrections (but ignoring 

resonance as the neutron star inspirals, the narrow set of points 
that fall within the resonance window Sa/R ~ 0.5(a/10i?)~^/^ 
should be ignored (see Lai 1994 and § 3.1). 
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Figure 7. Fractional residual (red solid line) with finite frequency corrections due to the non-static linear tide, as a function of 

orbital separation a/R for the neutron star model. The blue dashed line shows i.R)/R (s® plotted in Fig. 1) and indicates the location 
of linear resonances. The g-mode pair is f = 2, n = 1000 and (mi, m 2 ) = (2, —2), (2, 0), (2, 2) clockwise starting from the top left panel. 


linear damping) is then 


(WAB, VZH), we then have 


— (s — 


9j 






9j9^ 


-E 


(s — mp.Vt)^ — uj'^. 


(107) 


(r -b af - [r^ - uj^] - e^r 

+ ^lT^9g 




(108) 


where r = s — mpQ and a = {rup — mg)Q. There are two 
stable high frequency solutions near r ~ zLujp and two 
low frequency solutions near 


As in § 2.2, we can obtain an accurate estimate of the 
stability of this potentially large system of modes by con¬ 
sidering a two mode system consisting of a single p-g 
pair (we have verified this by numerically solving the 
eigenvalue problem given by equation 26 for large sets of 
strongly coupled modes; see § 5.4). Since K^pg ~ eupfuig 


r^-a±^u:l{l-ngg)-e^a^ (109) 

(cf. Equation 39). In the absence of finite frequency 
corrections {TZgg ~ 0) we recover the unstable solution 
of the incompressible limit in which the growth rate is 
r « 2ell if ujg < 2eVL (see VZH and § 2.2). If, however. 
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Figure 8. Magnitude of the residual Rgg as a function of g-mode 
frequency ujg/uo at a/R = 12 (/gw — 100 Hz; this is not near a 
linear resonance, see Fig. 1). We show results for g-mode pairs 
with £ = 2,4, 6 and mi + m 2 = 0. We find \Rgg\ oc e^{LOg/uJo)‘^, as 
shown by the magenta line (different values of mi + m 2 show the 
same scaling). 

the finite frequency corrections are significant and TZgg > 
1 , there is an unstable solution with a growth rate 

r « Wg yjTi-gg ■ (110) 

In Figure 8 we show \'R-gg\ as a function of ujg at /g^ — 
100 Hz (this is not near a linear resonance). We find that 
Tl-gg is well-fit by a relation of the form 

|R„I-av(^)’, (111) 

where the factor A = A(a;), which we describe in detail in 
§ 5.3, is a function of orbital separation and depends on 
and The w~'^ scaling is a consequence of the 
fact that the 0{ujg'^) and cancel in Equation (99) 

but the terms do not. The instability criterion 

TZgg > 1 is then equivalent to iOg < Aewo and the growth 
rate can be expressed as 

r « AetLiQ. (112) 

The lack of cancellation due to finite frequency correc¬ 
tions therefore leads to growth rates that are potentially 
a factor of « Awq/2 H larger than the incompressible limit 
growth rate. 

The sign of Rgg depends on the quantum numbers of 
the particular g-mode pair and in practice we find that 
there is no preference for either sign. In our numerical 
experiments with more realistic multi-mode systems, we 
find that the residual drives an instability whenever there 
are mode pairs with \Rgg\ > I. 

5 . 3 . Value of X(uj) 

In Figure 9 we show the magnitude of the residual fac¬ 
tor \{uj) as a function of a/R. The red lines show A 


for mode pairs with ^ = 2 ,n = 1000 , and {mi,m 2 ) = 

(2,-2), (2,0) and (2,2). The blue curve is Xh\^)^ 
same curve as in Figure I, and illustrates the variation 
of with a/R and the location of linear resonances. 
We find A ~ 0.1 — I in the regions outside the narrow 
resonance windows (as explained in footnote 2 , the re¬ 
gions near resonances where A has large peaks should be 
ignored). 

Although the fractional residual agg tends to decrease 
with increasing a/R (Fig. 7), we see in Figure 9 that 
the total residual Rgg = 2 aggK^( 2 )gg oc A^ is relatively 
constant with a/R (modulo and the resonance-to- 
resonance fluctuations). This is because the inner turn¬ 
ing point of the non-static tide, which occurs at the ra¬ 
dius where N{r) ~ w, moves inward to smaller r as 
oj = 2V, = 27r/gw decreases. The amplitude of the 
high-order g-modes are larger at smaller r, and therefore 
K^( 2 }gg increases significantly as /gw decreases (WAQB 
discuss this effect in the context of three-mode coupling 
of 5 -modes to As a result, although globally ca/N 

decreases, the finite frequency corrections remain signif¬ 
icant in the region that dominates the coupling (where 
oj « N{r)) and the residual only slowly approaches the 
static tide limit Rgg —)• 0{e^) as /gw —t 0. This slow 
convergence of the non-static linear tide to the static tide 
solution in local regions where N'^ is small is also noted 
in Terquem et al. (1998). 

We illustrate this point further in Figures 10 and 11. In 
Figure 10 we show, as a function of r, the radial displace- 
ment Xr of the {i, m) = (4,4) harmonic of the nonlinear 
tide and the cumulative integral d fo'' a/R = 5 

and a/R = 34, corresponding to/gw — 380 Hz and 22 Hz 
(neither of which is near a linear resonance). We assume 
an {£, m, n) = ( 2 , — 2 , 1000 ) self-coupled 5 -mode and also 
plot the horizontal displacement gt- At /gw — 380 Hz, 

( 2 ) 

Xr has an effective wavelength « R and most of the 

C 2 ') 

contribution to comes from large radii where Xr ' 

peaks (the crust at r > 0.9i? does not contribute because 
we assume it is not buoyant and therefore the 5 -mode dis¬ 
placement rapidly vanishes there). At /gw — 22 Hz, Xr^'^ 
has a short wavelength (^ R) and most of the contri¬ 
bution to comes from the region near the inner 

turning point of x?'^ (located at r ~ 0 . 2 i?), where the 
density and displacements of all three modes are large. 

In Figure 11 we show 2 K^(2}gg as a function of /gw. For 
5 -mode pairs that couple to the m = 0 and m = 4 nonlin¬ 
ear tide we find K^( 2 )gg cx f~J‘ while for pairs that couple 
to the m = 2 nonlinear tide K^( 2 )gg cx /g^. The former 
scale more strongly because the m = 0 and m = 4 non¬ 
linear tide is driven by two instances of the non-static 

linear tide {x^L 2 ,XmL- 2 } and {XmL 2 .XmL 2 }> respec¬ 
tively. The finite frequency corrections to the linear tide 
therefore enter the calculation twice. The m = 2 nonlin¬ 
ear tide, by contrast, is only driven by a single instance 

of the non-static linear tide {x^= 2 tX^=o}- 

Figure 9 also shows A for 5 -mode pairs with rii 7 ^ 72 - 2 . 
As the wavelength of the oscillatory component of 
decreases with decreasing /gw, it begins to couple well to 
daughter pairs with |ui — 712 ! < u^( 2 ), where 71 ^( 2 ) is the 
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Figure 9. Magnitude of the residual factor A (red solid line) as a function of ajK. The blue dashed line shows plotted in 

Fig. 1) and indicates the location of linear resonances. The p-mode pair is £ = 2, n = 1000 and (mi, m 2 ) = (2, —2), (2, 0), (2, 2) clockwise 
starting from the top left panel (the same pairs as in Figure 7). The points show A at odd integer multiples of ajK for £ = 2 pairs with 
m = 722 = 1000 (black points) and ni 7^ 712 with ni^2 = 1000 di 2 (magenta points). 


effective radial order of (i.e., the number of nodes). 
This is a consequence of linear momentum conservation 
and is discussed in WAQB in the context of three-mode 
coupling to the linear tide. At a/R < 5, n^( 2 ) ~ 0 and 
the coupling is much stronger for daughter pairs with 
ni = n 2 (black points). At larger a/R, n^( 2 ) > 0 and the 
coupling is also strong for pairs with ni ^ n 2 (magenta 
points). As a result, the effective residual Rgig 2 (and thus 
A) in the characteristic Equation (107) is, after summing 
over the « ^^( 2 ) strongly coupled daughters, larger than 
the ni = n 2 case by a factor of « n^( 2 ). This collec¬ 
tive driving of strongly coupled, multi-mode systems is 
described in WAQB and results in growth rates that are 
« .yn^( 2 ) faster than the single daughter pair estimate. 


5.4. Growth rate and the number of energy e-foldings 

In Eigure 12 we plot the growth rate T as a function of 
/gw, found by numerically solving the eigenvalue problem 
defined by Equation (26). We show results for (;-modes 
with £ = 2 and a range of parameter values (An, Ug, A) 
defined such that g-mode pairs in the range Ug ± An 
couple to each other with a residual factor A. Based 
on Eigure 9, A ~ 0.1 — 1 when the system is not near 
a linear resonance. Eor A = 0 we recover the incom¬ 
pressible limit result that T ~ 2en when ujg < 2ei2. Eor 
A ^ 0 and An = 0, we find growth rates T ~ Acwq 
when ujg < Aewo, consistent with Equation (109). And 
for the reasons given in § 5.3, we find that the effective A 
increases by a factor of « VSAn for An 0. In compar¬ 
ison to the A = 0 results, the finite frequency corrections 
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The number of mode energy e-foldings is approxi¬ 
mately given by 

rt r/gw 

ln{E/E,) ~ 2 / Tdt = 2 Tf-^df, (113) 

where /gw,i is the gravitational wave frequency when the 
instability first turns on and f /f = 2ta/^ (Equation 47) 
for the gravitational inspiral of two point masses. 

If we ignore finite frequency corrections (Tigg — 0) and 
damping, then T « 2eO (Equation 40) and 


HE/E,) = 


■fl/3 

^100 


fl/3 ■ 
^ 2,100 


(114) 


for M' = M and a neutron star with dynamical fre¬ 
quency Wo = lO"* rad s“^ (where /loo = /gw/100 Hz and 
All .2 = M./1.2Mq). In this case, the instability crite¬ 
rion is (jjg ^ 2efl, and the modes first become unstable 
at a frequency 

/ \ 1/3 


Figure 10. Cumulative coupling integral ggi^ '^) lines; 

divided by 10^e^) and the radial displacement of (i^m) = 
(4,4) harmonic of the nonlinear tide (blue lines; divided by R) 
as a function of radius r/R. The dashed lines are for a/R = t 
(/gw — 380 Hz) and the solid lines are for a/R = 34 (/gw — 22 Hz), 
neither of which are near a linear resonance. The p-mode in the 
calculation is a self-coupled ^-mode with (£, m,n) = (2,-2,1000) 
whose frequency LOg ~ 10~'^a;o. The dotted black line shows its 
horizontal displacement (divided by 10^-R). 


If ojg = IO“‘^wo (corresponding to 5 -modes with £ = 2, 
n ~ 1000) then by /gw = 1000 Hz the number of e- 
foldings is ln(if/£'i) ~ 11 , i.e., the modes grow by a 
factor of ^ 10 ^ in energy. 

If we account for finite frequency corrections, T « Aewo 
(Equation 112) and for M' = M, 


HE/E,) ^ AAM/T (^) 


■ -2/3 


f- 2/3 

JlOO 


(116) 



Figure 11. Magnitude of 2 K^(2)gg = E-ggjagg (divided by e^) 
as a function of /gw for g-mode pairs with £ = 2, n = 1000 and 
mi -I-m 2 = {0,2,4} (red, blue, and black solid lines, respectively). 
The black dashed lines show and scalings. 

yield significantly larger growth rates and the instability 
is triggered much earlier in the inspiral. 


In this case, the instability criterion is approximately 
ujg < Aewo, and the modes first become unstable at a 
frequency 

If Wg = 10“"^wo and A = 0.5, then the modes grow by a 
factor of ^ 10^° in energy by /gw = 1000 Hz. If Wg = 
lO^'^wo and A = 1 (corresponding roughly to the effective 
lambda for An = 2; Figure 9), then the modes grow by 
a factor of ~ 10®° in energy by /gw = 1000 Hz. 

In Figure 13 we show the growth in mode energy E 
as a function of /gw for a range of parameter values 
(An,ng,A). We find that depending on the values of 
(An,ng,A), the modes grow by factors of ^ 10^° to 
3> 10®° before the merger. Moreover, because the growth 
rates can be large relatively early in the inspiral, the 
modes can reach large energies already by /gw — 100 Hz. 

The linear damping rates of 5 -modes increase with de¬ 
creasing Wg and 7 g « ojg for Wg « 10“"*^a;o (§ 6.3). We 
find, however, that even allowing for a linear damping 
rate 7 g ~ Wg does not significantly limit the growth (see 
solid and dashed black lines in Fig. 13). 

5.4.1. Initial mode energy Ei 

What is the initial mode energy Ei and how many e- 
foldings are needed to reach dynamically significant en¬ 
ergies? A 5 -mode with i = 2 and m = 0 is linearly driven 
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Figure 12. Nonlinear growth rate T (in units of eo^o) of £ = 2 
p-modes (and the p-modes to which they couple) as a function of 
/gw for different values of the parameters (An, rig, X). 



Figure 13. Mode energy growth relative to an initial mode en¬ 
ergy Ei as a function of /gw for £ = 2 p-modes (and the p-modes 
to which they couple). The curves are for different values of the 
parameters {An,ng, X). The solid and dashed and black lines are 
equivalent except that the latter assumes that the p-mode has a 
linear damping rate 7 ^ = 


by the static linear tide to an energy 


« {elgfEo « 10 


-24 


10 


-4, 


OJq 




gw 


100 Hz 


Ei 


0 ) 


(118) 

where the linear overlap integral Ig ~ 0.3(wg/wo)^ 
(WAB). A low frequency (?-mode (wg <C w) with (£, m) = 


(2, ±2) is linearly driven by the non-static linear tide 
to an energy that is smaller than this by a factor of 
^ (wg/w)"^, where w = 2H is the tidal frequency. This 
yields Ei^g « 10“^®i?o for the values in the equation 
above. 

The near cancellation of three- and four-wave cou¬ 
plings Kj^gg -\- y) K^pgK^pg implics that in order to pre¬ 
vent the runaway growth of the (/-mode due to its self¬ 
coupling term K 4 gg, the many p-modes to which the lin¬ 
early driven p-mode couples must grow along with it and 
maintain an energy 

EpC^iKspgfEg^i^e^^ Eg. (119) 

For strongly coupled p-g pairs, kp ~ kg which for our 
neutron star model implies 


AgNCs 


cs --lO^Ag ( ) Wo. (120) 


-1 


Therefore for £ = 2 p-modes 


I 10-4^0 


-4 




gw 


100 Hz 


Eg, (121) 


i.e., the p-modes must maintain a much larger energy 
than the p-mode in order to preserve the three- and four¬ 
mode cancellation throughout the linear driving of the 
1 = 2 p-mode. Thus, by Equation (118) 


These estimates suggest that there are p-g unstable 
modes with initial energies Ei ^ 10“^° — 10“^°i?o at 
/gw = 100 Hz. For such modes, growth factors of 
E/Ei ^ 10^° are significant as they correspond to mode 
energies that are a substantial fraction of Eq, the star’s 
binding energy. Given that a p-mode with frequency 
ojg « IO^'^wq breaks at A ~ 10“^°Eo (i.e., |pgfcgPr| ^ 1; 
WAB), even growth factors ^ 10^° may be significant. 

We conclude, therefore, that the finite frequency cor¬ 
rections may lead to substantial mode growth prior to 
merger. Since each p-mode couples to many p-modes, 
there can be a very large number of excited modes 
10^). Determining what affect this may have on the 
inspiral requires an understanding of the saturation of 
the instability and is left to future work. 


6 . RESULTS WITH LINEAR DAMPING 

In this section we determine how linear damping influ¬ 
ences the stability of the tide to p-g coupling. In order to 
separate the effects of damping from the effects of finite 
frequency corrections, we will ignore the latter and as¬ 
sume that the cancellation given by Equation (36) holds. 
The characteristic Equation (35) is then 


- (s - uig. H) ^ -I- i^g. (s - rrig. H) -I- 


CO 


9j 


SpiQi 


Qi Pi 


(s — mp.fl) — i^p. (s — nip-fl) 
(s — mp.fl)^ — i^p. (s — mp^VL) — u 


9g, 


(123) 
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Figure 14. Instability growth rate P (in units of 2eQ) including 
linear mode damping as a function of gravitational wave frequency 
/gw for an -^ = 2 p-g pair. The mode frequencies are (ujp^ujg) ~ 
(2 X 10^, 10~^)c<Jo- The two upper lines are for p-mode damping 
rates of 7 p = 5 x 10^ s~^ and 10^ s~^. Both also assume a p-mode 
damping rate ')g = Ug. The lower line assumes no damping of 
either mode. The dotted black lines show y/'^p/AQ.. 


where we neglect 0 {e^) terms since they do not alter 
the stability. As in § 2.2, we can estimate the stability 
of this potentially large system of modes by consider¬ 
ing a simple two mode system consisting of a single p-g 
pair. Since K^pg ~ eujpjtOg, if we let r = s — mpft and 
a = {rUp — rUg) H, the characteristic equation becomes 
(cf. Equation 38) 

(r -I- a)^ — ijg {r a) — Ug 
- [r-^ - Z7pr] ~ 0. (124) 

We consider the solution to this equation for various 
limits of a and 7p. In § 6.1, we consider the case of the 
static tide (a = 0) with damping and show that all so¬ 
lutions are stable. In § 6.2 we consider the case of the 
non-static tide with damping. We show that p-modes 
that are strongly damped with 7p ^ H are unstable and 
grow, along with the g-modes to which they couple, at a 
rate T ^ eH (yp/H) . In § 6.3 we estimate the damp¬ 
ing rate of high-order p- and g-modes and argue that 
7p 3> H if the p-modes are above the acoustic cutoff fre¬ 
quency of the neutron star. In § 6.4 we show that for 
p-g pairs containing such strongly damped p-modes, the 
modes can reach significant amplitudes before the neu¬ 
tron star merges. 

6 .1. Static tide with linear damping 

By angular momentum conservation (i.e., the coupling 
coefficient selection rules Equations 17-19) only p-g pairs 
that satisfy |mp — mg\ = m couple to harmonic m of the 
i = 2 linear tide. Therefore for the case of the static 
tide TO = a = 0, Equation (124) gives, upon substituting 


[r^ - i-fpr - uil] 


r = —1/3, 

+ /3^[7p + Ig] + /3^[wp(l -I- e^) + Iplg] 

-l-/3[Wp(e^7p -I- 7p) -I- ypWg] -I- WpW^ = 0. (125) 

The Routh-Hurwitz criterion (Gradshteyn et al. 2007) 
for a quartic equation with real coefficients a 4 x'^-\-a 3 X^-\- 
a 2 x'^-\-aiX-\-ao = 0 shows that the solution is stable since 
Oi > 0 for all i, 0302 > 0401, and 030201 < 0401 -I- O3O0. 


6.2. Non-static tide with linear damping 
For the case of the non-static tide to = ±2, the low 


frequency solution to Equation (124) is 

+ 7g) 

(e^7p + 7g)^ 


=,-!> + -( 


± 


-g^- 


- e a[a + i{-ip- 7^)] 


1/2 


.(126) 


If the last term under the square root dominates, i.e. 

a;g,e^7p,7g, and 7 p>|a|, 7g, (127) 

then r ~ —a ± e(l -I- i)^J~orp^j2 and there is an unstable 
solution with a growth rate 

r « (128) 

Thus, the solution with damping is larger by a factor of 
« (7p/4r2)^/^ compared to the solution without damping 
or finite frequency corrections T « 2er2 (see Equation 40). 

To understand why large linear damping can result in 
a viscous-type instability that enhances the growth rates, 
note that the cancellation between three- and four-wave 
interactions requires that the phases of the p- and (/-mode 
oscillations be in strict relation. This relation can only 
be perfectly satished for the static tide (which is why 
the static tide is always stable). For the non-static tide 
without linear damping (nor finite frequency corrections), 
there is a phase offset in the equations of motion due to 
the Coriolis term ~ As a result, the growth rates 
scale as ^ w, as VZH showed (see also Equation 40). For 
the non-static tide with linear damping, the equations 
of motion have an additional viscous term ^ yw and, 
in analogy with the effect of the Coriolis term, we find 
that for large enough yp it can result in a growth rate 
~ (ypw)^/^ (Equation 128). 


6.3. Damping rates 


First consider the p-mode damping rate. By Equation 
(120), a (/-mode with ojg ~ 10“^wo couples to a p-mode 
with Wp ^ lO^AgWo- Such a p-mode is well above the 
acoustic cutoff frequency of a neutron star (wac ~ 10 ^Wq! 
WAB). Because it does not reflect at the stellar surface, 
it escapes in one group travel time across the star and in 
a region between r and r-1-Ar it has an effective damping 
rate 


27r 

tp(Ar)’ 


(129) 


where its group travel time across Ar is 



(130) 
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The size of the p-g coupling region determines Ar. If 
the driving rate is slower than the g-mode’s group travel 
time across this region, the driving is global and Ar ss R 
(WAB). In that case tp « R/cg « since Cg — 

1.5RuJo for a neutron star core. We therefore hnd a p- 
mode damping rate 

7p«I0^s"^ (131) 

This estimate motivates the value for jp used in the e- 
folding calculation below. In Appendix E we show that it 
maybe more appropriate to treat the g-mode driving as 
local rather than global. In that case, the coupling region 
Ar <C R and we estimate that the p-mode damping rate 
would be > 10® s“®. 

In WAB we found a p-mode damping rate 

« 3 X lO-^A^Tg-^ s-i, (132) 

where Tg is the core temperature in units of 10® AT. The 
minimum p-mode frequency such that ujg > 7^ is there¬ 
fore 

Wg « O.TA^/^Tg”^/^ rad s"^ (133) 

In practice, we find that the modes are p-g unstable even 
if LUg < 7 g, although in this regime the g-modes are more 
correctly treated as traveling waves rather than global 
standing waves. In any case, this motivates our choice of 
ojg = lO^^wo in the estimates below. 

The damping rate calculations above assume that the 
modes are oscillating at their natural frequencies. How¬ 
ever, p-g pairs will oscillate far from their natural fre¬ 
quencies when they are unstable to nonlinear p-g driving 
by the tide (see Equation (126)). The p-mode, in partic¬ 
ular, is forced to oscillate at a frequency ^ ujp. Although 
not entirely clear, we might still expect a p-mode whose 
natural frequency is above the acoustic cutoff to have a 
large damping rate that is given approximately by Equa¬ 
tion (129). This is because once the p-mode leaves the 
coupling region in the core, it is no longer unstable and 
forced to oscillate at the nonlinear driving frequency. In¬ 
stead, upon entering the crust, it will begin to oscillate at 
its natural frequency. It therefore does not reflect at the 
stellar surface, and thus its effective damping rate is ap¬ 
proximately the inverse sound crossing time. A caveat is 
that its group velocity within the coupling region might 
be less than Cg. Therefore, perhaps only the outgoing 
acoustic waves excited near the outer coupling region 
have damping rates as large as the estimates above. 

6.4. Number of e-foldings with damping 

Eigure 14 shows the p-g growth rate accounting for 
rapid mode damping. In order to obtain these results, we 
numerically solve the eigenvalue problem (Equation 26). 
We assume that the p-g pair has frequency (u!p,ujg) ~ 
(2 X 10®, 10“‘*)wo- Based on the estimates in § 6.3, this 
suggests p-mode damping rates 7p « 10® s“® and p-mode 
damping rates 7g ss ujg. We find that the numerical 
calculation agrees well with the analytic estimate given 
by Equation (128). The growth rates are significantly 
larger when damping is included and the instability is 
triggered much earlier in the inspiral. 


By Equation (113), the growth rate of Equation (128) 
results in an energy e-folding 


HE/E,) cs l98MfT 



P-1/6 

Ji,100 


P-1/6 

/lOO 


(134) 

for M' = M and ujq = 10^ rad s The instability cri¬ 
terion is ujg < €^y2nH, and the modes become unstable 
at a frequency 



(135) 


If the instability begins at /gw,i = 60 Hz then the mode 
energy grows by a factor of E/Ei ~ 10®® by /gw = 
1000 Hz. 

We show in Appendix E that the p-mode driving might 
be in the local regime, in which case 7p 10® s“®. The 
modes would then grow by E/Ei ^ 10®®. 

7. CONCLUSIONS 


We analyze the stability of the tide in coalescing binary 
neutron stars to p-g mode coupling. Previous studies did 
not account for either the relevant four-wave interactions 
(WAB) or the compressibility of the linear tide (VZH). 
In order to account for both effects, we first solve for 
the non-static linear tide and thereby obtain the 

finite frequency corrections to the static tide (such as 
compressibility). We then use this solution to evaluate 
the three- and four-mode coupling coefficients. We find 
that the finite frequency corrections are significant and 
undo the near-exact cancellation between the coefficients 
found in the incompressible limit. As a result, the insta¬ 
bility begins earlier in the inspiral and the p-g growth 
rates are ~ Wq/w larger than those of the incompressible 
limit. We find that the unstable p-g modes can poten¬ 
tially reach significant energies well before the neutron 
star merges. 

In a separate analysis, we investigate the effects of lin¬ 
ear damping on the p-g instability (ignoring finite fre¬ 
quency corrections in order to disentangle the two ef¬ 
fects). For small linear damping rates, the damping acts 
to slow the mode growth. However, if the p-mode damp¬ 
ing rate is sufficiently large, it induces a viscous-type in¬ 
stability and the p-g growth rate increases with increas¬ 
ing 7p. We find that p-modes above the acoustic cut¬ 
off frequency of the neutron star have effective damping 
rates 7p > ujq and increase the instability growth rates 
relative to the inviscid limit by a factor of ~ (wo/w)®/^. 
This viscous mechanism can therefore also drive modes 
to significant energy before the merger. 

In this study we focused on the stability and growth 
rate of p-g coupling. We did not attempt to solve for the 
instability’s saturation and therefore did not determine 
its influence on the tidal dissipation rate. Thus, we do 
not know to what extent the instability affects the phase 
evolution of gravitational waves from coalescing neutron 
star binaries. Nonetheless, its early onset (/gw ~ 50 Hz), 
rapid growth rates, and large number of excited modes 
H 10^): suggests the instability’s impact could be signif¬ 
icant and motivates the study of its nonlinear saturation. 
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APPENDIX 

A. STABILITY OF A STATIC TIDAL FIELD 


Consider a neutron star that is deformed by a static tidal field. We show here that a short wavelength perturbation 
to this static, tidally deformed neutron star is stable to all orders in the tidal factor e = {M'/M){R/a)^ as long as 
e 1 (we thank G. Ogilvie for bringing this argument to our attention). 

The response of a static, tidally deformed star to a small perturbation ^ is given by 


52 ^ 


-V6p - 6pV4) - pVS<l>, 


(Al) 


where 5p, 5p, and 54> are Eulerian perturbations and the gravitational potential </> is that of the star plus the (time 
independent) external tidal potential. Since the background is taken to be the tidally deformed star (not a spherically 
symmetric star), the background quantities p,p, etc. account for the exact static tide response to all orders in the 
tidal factor e. If we assume that p and p vanish at the surface of the star and the perturbation varies with time as 
^ ~ e*®*, the energy equation of the perturbation is 


s 


1 



I 

4^ 



\V6(l)\^dV 



+ it ■ Vp) (4 • A) 


dV, 


(A2) 


where V is the volume of the tidally distorted star and A = V Inp — (I/ri)V Inp is the Schwarzschild discriminant. 
Since the perturbation is assumed to be small, we only account for its energy to lowest (i.e., quadratic) order. 

The three contributions to are: the perturbed gravity term (negative and therefore destabilizing), the perturbed 
pressure term (positive and therefore stabilizing), and the buoyancy term, whose influence on the stability we now 
consider. In a spherical star, (|* • Vp)(^ • VA) = pN‘^\^r\'^, and the buoyancy term is stabilizing if > 0, where 
= Vp ■ Afp is the Brunt-Vaisala buoyancy frequency. Consider now a star deformed by a static tide. Since the 
hydrostatic equilibria are barotropic, the vectors Vp and A are parallel and (|* • Vp)(^ • VA) = pN^\^g\'^, where 
now pN"^ is that of the tidally deformed star and is the component of ^ along g = — V((). If we start with a stably 
stratified spherical star and apply a static tidal deformation, we still have > 0 as long as e < I. This is because the 
Eulerian perturbation to due to the static tide is SN^/N^ ~ e (see, e.g., Burkart et al. 2013). The buoyancy term 
in the tidally deformed star is therefore stabilizing if the spherical star is stably stratified, as is the case for neutron 
stars. 

For short wavelength perturbations, such as those we consider in the main text, the perturbed gravity term 6 (f> in the 
energy equation should be unimportant (i.e., it is appropriate to work in the Cowling approximation). We therefore 
conclude that a short wavelength perturbation to a static, tidally deformed neutron star is stable to all orders in e. 


B. USEEUL EORM FOR THE SUM RULE INTEGRAL 


We derive here an expression of the sum rule integral (Equation 87) that makes its similarities with K 4 cd more 
apparent. Assuming the Cowling approximation, the linear operator acting on some vector ^ can be written as 
(Reisenegger 1994) 


ZiM=_v 

p 




{V-i)g 


V[clV-i + g-^] + 



(v-Os 


(Bl) 


(the two restoring forces are the pressure variation 6 p and the buoyancy associated with N'^). For any two vectors ^ 
and g, we want to evaluate 


fj^[n] = / d^xp^- (v [clV -g + g-g] H- ■g)g 


(B2) 


We can use integration by parts to simplify the first term 


J d^xpi ■ V [c^V ■g + gg\= j dS ■ pi [c^V g + gg]- j d^xV ■ {pi) [c^V g + gg] 

= - J d^x {pcliV ■ i) {V ■ g) + p{V ■ i)g ■ g + i ■ Vp[clV ■ g + g ■ g]} , (B3) 

where the surface term vanishes because we are interested in vectors that satisfy the surface boundary condition 

{ag ■ r + C2V • g)^^ii = 0, (B4) 

where ci and C2 are constants (Equation 5 in Reisenegger 1994). We then have 

J d^xi- f4[g] = - J jrip (V • ^) (V • 77) + P (V • ^) g • J7 -1 • (p-^ (V • rj) g - Vp [c^ V ■ g + g ■ g]^^ .(B 5 ) 
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Using the definition of , the integrand term within large parentheses simplifies to 

iV^ dlnp] 

9 dr 


pclV ■ T] 


dp dp 

9^Vr=Pg" -V+g-Tdr 
dr dr 


and we find 


J d^x$, ■ /i [rj] = J d^ 


dp, 


-Tip (V • 0 (V • 77) + pg • rj + 77^V • |) + 

Note that this is symmetric in ^ -f-)- rj, as expected given that fi is an Hermitian operator, i.e., 

J d^x$, ■ fi [rj] = J d^xfi [^] • 77. 


(B6) 

(B7) 

(B8) 


Equation (87) can therefore be written as 


E —J- / <i’ 

&e{p} 


1 - 


2w^ 


tP*ac • fl bPadl - 7^ • fl [Cad] + ^’ad ' fl [Cac]) 


1 

Eq 


J d^X -Tip (V • rae) (V • 7^,,) + pg {rac,r^ ' ad + C’ad,.V • ^ac) + 9^rac,ri’ad,r 
-^ [ d^X^ [2a;^ (V • ^ac) (V • ^Pad) + ^0 {(V • ^ac) (V • Cad) + (V • CL) (V • 7^,,)}] 

^0 J \‘^b) 

I d^X -Tip (V • 7^:J (V • 7^,,) + pg {^acr^ ■ i’ad + P’ad,rV ■ ^ac) + g^rae,ri’ad,r 

^^ac,r dCad,r ac,r dipad,, 


1 

Eq 


riP 

K) I 


2 ^^ac,r d'lljad,r 2 

-dkr 


dr dr 


dr dr 


(B9) 


where we dropped terms smaller than given that (w^) oc uj,, ~ and V • ~ dtpac,r/dr 


0{uj,, ^). In terms of the 0{uj,, function 


Zac{r,0,(p) = r 
we have, after integration by parts. 


V • 7/^ac - 


dlpac,r 

dr 


— ‘^P^ac^r E 


sin 9 


d {sin Oipacfi) dipac, 4 > 
dO ^ d(p 


(BIO) 


f ,3 f j3 f d{tpac,r'kad,r) 4’ac,r^a,d + 4’ad,rZa 

j d-^xpg ■ 7/>^^ + 7/-ad.rV • ^ac) = / d^^xpg I - ^^^ - 


We therefore have our final result (Equation 88) 
KstcK^bd-^ ( d^a;|rip[ 


dg dp 2pg 


'4’ac,r'4’ad,r 


‘Pac,rZ'ad + '4’ad,rZa 


(Bll) 


be {p} 


E, 


0 J 

PI 

r 

Tip 

{^D 


(oy-4) (4^^) + {d^^-ks) + (c*4d) 


'^ac.r^ad H“ '4^ad,r ^ac ( 2 -|- 


ding\ 


dim 


'4^ac,r'4^ad,^ 


2 ^^ac,r fi7/>ad ,r . 2 f ^^ac,r d(^ad,r , ^Cac.r ^V^acZ,r 

2CJ ---^ ^0 


\ dr dr dr dr 


+ e 0{uj^ ), 


(B12) 

(B13) 


dr dr 

where we wrote V • in component form 

V • tA,, = V • [(c . V) a] = = a^c7 + 

and used the fact that a = 

C. FOUR-MODE COUPLING COEFFICIENT K^bcd 

Here we give the form of the four-mode coupling coefficient Kabcd that we use for numerical calculation. We derive 
it using the same general approach taken by Wu & Goldreich (2001) and WAQB to derive a numerically useful form 
for the three-mode coupling coefficient Kabc- We find that the coefficient can be broken into five pieces 


^abcd — 


/ + // + /// + /y + y 
^0 ’ 


(Cl) 
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and can be written in terms of angular integrals, ^r, ^h, V • and 


A U — ^(1) r(l) /,(2) f(2) l(3) A3) 

“ '‘‘ab Jab + 'Gfa Jab + '''ab Jab > 


( 02 ) 


where 






, A h 

- rx n btr-Oh ^ Cihar q Ctra', 

or or r^ r^ r 


2 2 ' ' 2 
Or dbh , br dtth 1 d{ahbh) 


r dr r dr 


dr 


YaYb, 

VYa • vn, 




as follows: 

-f = Tabcd [ drr'^pGiV ■ aV • 6V • cV • d 


11= drr^pG2 


V aV b 


3 3 


y'-cd 


I /,(2) p(2) I 1,(3) p(3) 

bed -I- il-cd ^ab,cd ""cd ^ ab,cd 


{permutations of abed} 


III = / drr' 


'Yip EE| ^ab ^cd ^ab,cd + /^S4dS£L + /^id^bf^^2L 


i=i ]=i 

IV = I drVip-l — 4r^V • aV • 5V • cV • dT^ 


abed 


+r^V • a 


V-6 


_ dbr 

+V-a( — 

or 


{hi^X,,d + + {6 o c o 4 

(2Crdr Crdh-^d dr-Ch-^c) ^abed “t” ^hdh {^c^d^abed F, 


(3) 

ab.cd 


F {Cr C/i) (^dfiSa^bc,d 44if/,fcc^ ^'>^^ac,bd 

+{permutations of abcd}|’ 

T/_ [ i h rl ^ h rl ( 1 ^ {^(22) , p(22) , p(22) 

^ ~ J '^^P^rbrCrdrr Tabcd + ClhbhChdh ( r dr ) V ^acM + ^ad,b. 

+ ~ r ^®'’&rCbdbF’^j,2 + ‘^rbhCrdhPj^Ji^^ + arbhCudrF^J^^^ 


F c (i} 


+ahbrCrdhF^J^^j^ + ^hbrCfidrFl^J^^^ + ahbhCrdrFl:J^^^ 


(03) 

(04) 

(05) 

(06) 

(07) 

(08) 


(09) 


(CIO) 


Here the divergences are with respect to the radial functions. The angular integrals are defined as 


Tabcd = J dnYaYbY.Yd, (011) 

Fald^ J dllY^YFI, (C12) 

FSed^ J d^Ffcd, (C13) 

Sa,bc,d = J dOY,V,nV^ Y,V, VWd + AFFc- (C14) 

Paired subscripts not separated by a comma are symmetric in those indices. Note that = Tabcd- The integrals 

(‘2'\ 

Tabcd, Tab cd' Tab cd' Sa,bc,d are, respectively, analogs of the three-mode integrals Tabc, Tt,, 14, and S defined in 
Wu & Goldreich (2001); the integrands of the former equal Ya times the integrands of the latter. We compute these 
integrals numerically in our calculations. 
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( 22 ) 

C.l. A relation between Sa,bc,d and E^a,b cd 

The differential expression under the integral in the definition of Sa,bc,d can be written as 

v.nvw.v.vWd = V, (vWdV.n) - v.idVWeVjV^n. 

Therefore 

SaM,d + Sa,cd,b = Jdn YaV, ,Yb) VWe + 

= I dnv^Y^v^Yb [KlYaX, - + a^f^J, 

_ A 2 77 ( 2 ) _i_ A 2 77’(2) I A ^ _ Z7’C22) 

^^b"ab,cd ' ac,bd ' ^^d^ad,bc ^ad,bc'> 

where we used integration by parts and the relation V^F = —h?Y. We thus have 

p( 22 ) o _ ^( 22 ) Q 

^ac,bd CJojM.c — njad,bc 'Ja,bc,d- 


(C15) 


(C16) 

(C17) 


D. EXPRESSIONS FOR THE 0{w-^) CONTRIBUTIONS TO Kigg + E^Sps-fCspg 

In this appendix we provide expressions for the 0{uj~‘^) contributions to K^gg + ^ K^pgR^pg. In § D.l we consider 
3K^(i)^(i)gg + ^3pg^3pg- We then describe two methods for calculating the nonlinear tide which we need in 

order to evaluate 2K^(2)gg: (§ D.2) as a sum over modes and (§ D.3) by directly integrating the inhomogeneous equation 
of motion. 

Once we have we can calculate K^( 2 )gg similar to how WAQB compute the three-mode coupling coefficient Kabc- 
Note, however, that we cannot use the final expression for Kabc in WAQB (A55-A62) since their analysis assumes that 
all three modes are non-radial modes {£ 7 ^ 0) whereas £ = 0 is one of the harmonics of Instead, we proceed 

as in VZH and use their Equations D10-D13. When calculating x*'^^ by direct integration (§ D.3) we must account 
for the inhomogeneous term in the equation of motion, which yields, in addition to VZH’s Equation DI0-DI3, an 
inhomogeneous coupling term 

j “ ^hdh) fr,y^(2)F^(2), (Dl) 

where c and d label the g-modes, 7 Q( 2 ) is the three-mode angular integral and fr,x<.^) is tbe radial component of the 
nonlinear tide driving force (given by Equation D34 below). 


D.l. 3K^(l)^(l)gg F F^vgF^na 

Starting from Equations ( 88 ) and (95), we decompose the terms that enter at 0{uj~'^) into eight separate pieces 


iKaacd + ^SbcK^bd = I + H + HI + IV + V + VI + V11 + VIII, 

b 

where a labels the linear tide (whose frequency is u;), c and d label the ^-modes, and 


1=^1 d^xT.p 

Fo 


3^t-. ^ I (cljal) (d’^a-^sk) + ) (df^affc) - V • a det | (a + c + d)^ 

I 


II = 


F, 


0 


d^xrip{c^aV) 


///= - 


^ f 1-^ P9 f , 

I ^ i H^ad.r^ac ( ^ dTn^ / 


En 


IV = - 




Fn 


K) I 


2 dl/’ac.r dV’ad.r , 2 / ac,r dpad,r , drj ac,r dtpad, 


2 ^^^ d 7 ' d 7 ' +‘^0 


dr dr 


+ 


dr dr 


V = -- 


2F( 


J d^xTip + 2^ ((V • af Fyb + 2 (V • a) (V • c) a(yb + 2 (V • a) (V • d) a^cf, 

d^xrip{algal) (cf.dffc) 


VI =- 


2Fn 


E//= ——b / d^a^Tip (V • c det (a + a + d)b 
2Fo J \ 


V • d det 


{a + a + c)b ^ 


(D2) 

(D3) 

(D4) 

(D5) 

(D 6 ) 

(D7) 

(D 8 ) 

(D9) 
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VIII = / (fxpala^c'^d’^(j).ijks (DIO) 

2Eo J 

In § 3.5 we show that term I is O{oj~'^) due to cancellations between its individual terms. From the definitions 

of and (Equations 71 and 72) and the four-mode angular integrals defined in Appendix C, we can reduce the 
portion of these terms to radial integrals. With the double sum over all harmonics of the linear tide mode a 
implicit, we find: 


J “ ^hSa,cd,a - + (^r - an) {e^^cM “ '^'a.ad.c 


d{V -a) 


{ar-ak)[hlFl^l, + klF^Xc 


+ripV • + A^A^T.,,, - ^ - S,,ac,d - 

^ ^ ^ (a^ [EiJ,, - - EiJ,J - 

+ ^ [2AlAlAlT,,,, + - A^ - A^ ) | c,4 (Dll) 


TTT--— ( ^^^ S°-r- 0 -h) (\, ding] ^( 22 ) _A 2 p( 2 ) _A 2 p( 2 ) 

TA I _ (i-r I 4 ^ ^ ^ac,ad ac,ad ^d^ad,a 


( fD 2) , p(32) _ L a2 _ dlng ^ (22) 

I -^ac,ad +'^ad,ac ^ac,a 


^ y dr(a, - a,,) (w" + co'^) « - a),) 

r dlnFi _ /ainFi\ 1 „7A2p(2) , A2p(2) \ 1 dcn ddh 

f ^ dlnp 51np jj ad,ac)\^ 

^ ^J 

-2p5 {ar - an) (aIF^^, + A^E®,,) ] + 4 [rFip + 2^ (V • a)^] e2,,| 


dlnFi _ plnFW ' 
dlnp \ 91np ) g 


2 V • a 
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D.2. Calculating with evaluated as a sum over modes 

Expanding the nonlinear tide as a sum over modes 

Equation (7) becomes 


Uc 


EEt^rxi^’+E KabcXb X. 
. b m be 


Xa^' + i^lXa^ = [Va{t) + Ka{t)]* = UjI 
where (cf., Equations 20-22 in WAQB) 

K(t) = J d^xpi^ ■ ■ V) 


(1)^(1) 




VC/, 


Ka{t) = ^ J /2 




(D19) 


(D20) 


(D21) 


(D22) 


The amplitudes of the nonlinear tide expansion are therefore given by the steady state solution to Equation D20 (cf., 
Equation 59 in WAQB) 

W>1, (K + Ka) i 






- {rria^)^ 


For the static tide (rua = nib = me = m = 0), 


r«=C/„, = - 


En 


d^xpC-^U 


and the solution is (cf. VZH Equation 23) 


X, 


i^) =Va+Ka = ^ (^abUb + E ^abcUbUej . 


(D23) 


(D24) 


(D25) 


In the more general case that includes the non-static linear tide (C = 2, m = {—2,0, 2}), we first write 

x^^\x,t) = X-2(x,t) +Xo(x,t) + X2(x,t) = e[VE2oXo(®) + 2IT22X2(®) cos(2nt)] • (D26) 

The harmonics that contribute to x*'^^ are £a = {0,2,4} and ma = {0,±2,±4}. The coupling coefficients for the 
different values of ma are then 


(Ea+iXa)„^^0 = e" [W"20f^20( 


^axo ^ ^^axo 


^ifk)%0.0]+2IT22fE22( 


7 ~ + 

aoX2 w ^'^aX2 


, CH) \ rr. 

+ '^aX2X2 ] ^[0.2-2] 


{Va+Ka)^^^^^ = e^W2oW22[ 

(K + = £"1^221^22 ( 


^axo ^ ^^axo 


J ~ 4- 2k^^^ 
^ax2 ^ ^^ax2 


^ax2 ^ ^^ax2 


+ 241^x2) ^[2.0.-2], 


, CH) \ rr. 

+ '^aX2X2] %-2-2], 


(D27) 

(D28) 

(D29) 


where T is the three-mode angular integral, which we label here by the values of m (the three C values are £a for the 

nonlinear tide and C = 2 for each linear tide). The coefficients Jab, K^Jb ^ and are as defined in Appendix A of 
WAQB (the hat symbol indicates that the angular integral T is factored out). 

D.3. Calculating K^( 2 }ed with x*'^^ determined directly from the eguation of motion 

Rather than expand x^^^ as a sum over modes as in § D.2, we find that it is numerically more accurate to instead 
directly solve the inhomogeneous equation of motion 




VC/. 


(D30) 


We do this by first solving Equation (6) for the linear tide X^^^- By angular momentum conservation (in the form of 
three-mode coupling between the linear and nonlinear tide), we have that the inhomogeneous driving term and x*'^^ 
both oscillate as where uj = m£l and m = (0, ±2, ±4} is a harmonic of the nonlinear tide. The time dependence 
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therefore cancels out and the equations reduce to a boundary value problem involving a pair of coupled linear ODEs 
in the radial direction with inhomogeneous driving terms. We find (cf., Pfahl et al. 2008) 
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dlnr 
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= - 3 h/i 


gklr 


gir 


c2 

dlnm 
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^ ■ dh^ + — ) 1/2 + - 


g/r) 
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pg 


where, letting a and b represent harmonics of and c represent a harmonic of 


Cr ^ I 2t7 

2/1 = —, 2/2 = —, tp = gcr - ■ c, 

r gr 


and the radial and horizontal driving forces are (see also Equations D39 and D40 below) 

fr = J2 [ dnY:{f^ [a,b]-p{a-V)VUt}-r, 


ab 


= dn ■ {/2 [a, 6] - p (a • V) . 


(D31) 

(D32) 

(D33) 

(D34) 

(D35) 


ab 


Here we use the notation Ut to indicate harmonic b of the tidal potential (and not the linear driving coefficient found 
in WAQB). Equations (D31, D32) assume that oj = WcH ^ 0- If a; = mcH = 0 (i.e., the static nonlinear tide) and 
£c 7 ^ 0, then the nonlinear tide is instead given by the solution to 


ijj 

dpj 

dr 


rfh 
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(D36) 

(D37) 


The first equation gives i/>(r) and we can use that result to get Cr from the second equation. If = 0 then Equation 
(D37) still holds and the equation for dcr/dr is, by the definition of the divergence. 


dcr 

dr 


— Cr + V • c = 
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Starting from Equation (137) in Schenk et al. 2002, the radial and horizontal driving forces can be written in terms 
of the linear tide displacement a and b as 
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(D39) 

-gJ 


(D40) 


At the center we impose the regularity conditions Cr = £cCh and at the surface we require the fluid to be hydrostatic 
by imposing Ap = 0. We solve the equations by shooting from these boundaries to a fitting point at r ~ i?/2 (see, 
e.g.. Press et al. 1992). 
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E. DAMPING RATE OF p-MODES IN THE CASE OF LOCAL DRIVING 

In § 6.3 we assume that the g-mode driving is global. However, if the driving is local then the coupling region 
Ar ^ R. In the calculation below, we show that in that case the p-mode damping rate is much larger than that in 
the global case (i.e., 7 ^ 10® s“^). The (/-mode driving is local if 

^tgiR) > 1, (El) 


where E ~ is the nonlinear growth rate and tg{R) ~ krR/ujg is the (/-mode’s group travel time across the star. 

For an equal wavelength p-g pair, kr — ujpics — Wp/l.SRwo, and the condition for local driving becomes 


- o X iU /loo' 


VlO® s- 


(E2) 


For ojp > Wac ~ lO^Wo, wavelength matching implies ujpjijjg > 10®and driving is local if /gw ^ 30 Hz. 

To estimate 7 ^ in the local driving regime, assume Ar is determined by the distance the (/-mode travels in a nonlinear 
growth time, i.e., equate the (/-mode’s group travel time across Ar to its inverse growth rate 


tg{Ar) 


-^-1 


Solving this equation for Ar assuming an equal wavelength p-g pair yields a damping rate 
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We thus see that 7 ^ ^ 10® s ^ in the local driving regime. 


(E3) 

(E4) 
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